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SUMMARY 

The  behaviour  of  a vortex  sheet  embedded  in  an  irrotational  steady  flow  and 
springing  from  a separation  line  on  a smooth  body  is  studied  in  the  neighbourhood 
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1 INTRODUCTION 

Classical  aerodynamics  is  based  on  regarding  the  flow  past  a streamlined 
body  at  a small  angle  of  incidence  as  inviscid  and  irrotational  outside  a thin 
region  which  comprises  the  boundary  layer,  surrounding  the  body,  and  the  wake, 
extending  downstream  of  it.  The  fluid  in  the  boundary  layer  and  wake  is  charac- 
terised by  its  acquisition  of  rotation  through  the  action  of  viscous  forces  and, 
at  high  Reynolds  numbers,  by  its  turbulence.  In  calculating  a first  approxima- 
tion to  the  external  flow,  displacement  effects  of  the  boundary  layer  and  wake 
are  neglected,  while  the  net  rotation  in  the  wake  of  a lifting  surface  is  repre- 
sented by  a planar  vortex  sheet  whose  strength  is  determined  by  the  Kutta- 
Joukowski  condition  of  finite  velocity  at  the  trailing  edge. 


When  the  angle  of  incidence  is  no  longer-  small,  the  boundary  layer  may 
leave  the  surface  of  the  body  well  upstream  of  the  trailing  edge.  This  is  known 
as  'separation',  although  the  process  is  not  essentially  different  from  the 
departure  of  the  boundary  layer  from  the  trailing  edge  in  the  classical  aero- 
dynamic model.  In  the  separated  flow  which  results  from  boundary  layer  separa- 
tion, rotation  is  present  not  only  in  the  boundary  layer  and  wake  downstream  of 
the  body  but  also  in  the  external  flow  which  surrounds  it.  Modelling  such  flows 
presents  difficulties  of  two  kinds:  practical,  computational  difficulties  of 
representing  the  rotation,  and  conceptual,  physical  difficulties  of  accounting 
for  its  origin.  In  general,  these  difficulties  are  compounded  by  a degree  of 
flow  unsteadiness  very  much  greater  than  that  of  a turbulent  boundary  layer 
upstream  of  separation.  This  unsteadiness  leads  the  aircraft  designer  to  avoid 
such  flows  as  far  as  possible,  so  there  is  some  justification  for  their  neglect 
by  the  aircraft  aerodynamicist.  However,  if  a separation  line  is  highly  swept, 
the  separated  flow  associated  with  it  is  frequently  steady,  particularly  when  the 
position  of  the  separation  line  is  fixed  by  a salient  edge  on  the  body.  The 
prime  example  is  the  flow  over  slender  wings  at  incidence,  with  the  formation  of 
coiled  vortex  sheets.  Such  a flow  is  eminently  applicable. 


Considerable  success  has  been  achieved  in  modelling  flows  with  highly-swept 
separation  lines  along  salient  edges  by  a straightforward  extension  of  the 
classical  aerodynamic  model.  The  rotation  is  still  represented  by  a vortex 
sheet,  the  displacement  effects  of  the  separated  flow  are  ignored,  and  a Kutta 
condition  is  applied  at  the  salient  edge;  on  the  other  hand,  the  location  and 
form  of  the  vortex  sheet  can  no  longer  be  assumed  a prioTi.  Most  work  so  far  his 
made  use  of  the  framework  of  s lender-body  theory  to  reduce  the  computational 
complexity,  but  it  has  recently  been  demonstrated1  that  the  approach  can  be 
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carried  through  more  generally*  in  particular  for  fully  three-dimensional 
incompressible  flow. 

Highly-swept  separation  lines  also  arise  on  smoothly  curved  bodies*  for 
instance  on  the  noses  of  aircraft  at  high  incidence*  on  guided  weapons  and  on 
up-swept  rear  fuselages  of  aircraft.  Parts  of  the  leading  edges  of  Concorde  are 
appreciably  rounded  and  the  leading  edge  of  the  US  Space  Shuttle  is  well  rounded. 
Elsewhere  on  slender  wings  highly-swept  separation  lines  are  found  on  the  upper 
surface  beneath  the  primary  vortices  and  may  arise  on  the  shoulder  of  a deflected 
leading-edge  flap. 

It  is  therefore  of  considerable  interest  to  attempt  to  extend  the  coiled 
vortex-sheet  model  of  separation  from  a salient  edge  to  model  steady  separation 
from  a smooth  surface.  If  separation  originates  at  a pointed  apex*  as  for 
instance  on  an  ogival  nose  at  high  incidence*  the  separated  flow  can  be  repre- 
sented by  a coiled  vortex  sheet  which  takes  a conical  form  at  the  apex;  but  if 
the  most  upstream  point  at  which  the  flow  separates  is  an  ordinary  point  on  the 
surface,  as  on  an  ogival  nose  at  lower  incidence*  its  initial  form  is  not  clear. 
Apart  from  this  initial  problem*  we  may  expect,  based  on  experience  in  calculat- 
ing flows  with  salient  edges,  that  the  boundary  conditions  on  the  vortex  sheet 
(that  it  is  a stream  surface  with  no  discontinuity  of  pressure  across  it)  will 
define  the  sheet  completely  once  the  separation  line  from  which  it  springs  has 
been  specified.  (There  may  be  more  than  one  solution  fcr  a given  body,  onset 
flow  and  separation  line,  as  found  in  Ref  2,  but  each  is  well-defined  by  the 
conditions.)  The  position  of  the  separation  line  must  be  supplied  to  the 
inviscid  model;  it  might  be  obtained  from  experiment  or  through  an  iteration 
between  the  calculation  of  the  inviscid  flow  and  the  calculation  of  the  boundary 
layer  on  the  body  upstream  of  separation.  The  background  to  these  introductory 
remarks  is  given  at  greater  length  in  a recent  review"^. 


It  is  in  the  immediate  neighbourhood  of  the  separation  line  itself  that 
the  main  differences  in  behaviour  are  likely  to  arise  between  separation  from  a 
salient  edge  and  separation  from  a smooth  surface.  The  present  analytical  study 
was  undertaken  in  order  to  shed  some  light  on  this  behaviour  and  to  provide  a 
guide  for  choosing  appropriate  local  representations  in  subsequent  numerical 
work.  The  results  obtained  define  the  possible  types  of  behaviour  of  the  vortex 
sheet  close  to  the  separation  line  and,  through  the  associated  pressure  fields* 
raise  questions  concerning  the  eventual  matching  of  the  inviscid  model  to  a 
boundary  layer  calculation. 


TWX  . 


In  section  2 it  is  shown,  with  some  minor  assumptions,  that  the  vortex 
sheet  leaves  the  surface  tangentially,  ie  that  at  almost  all  points  of  the 
separation  line  the  tangent  plane  to  the  sheet  is  also  a tangent  plane  to  the 
body.  The  argument  is  elementary,  but  is  is  set  down  to  provide  a basis  for 
what  follows  and  to  make  explicit  the  assumptions  involved. 

For  the  rest  of  the  study,  the  assumptions  of  slender-body  theory  are  made. 
This  provides  a manageable  analytic  framework  and  is  likely  to  be  appropriate  in 
most  situations  involving  highly-swept  separation  lines.  The  method  of  analysis 
employed  is  described  in  section  3 and  applied  in  detail  in  section  4 to  the 
problem  of  a conical  vortex  sheet  leaving  a semi-infinite  plane  wall.  It  is 
shown  that,  if  the  distance  of  the  sheet  from  the  wall  increases  like  y11  , 
where  y is  the  distance  along  the  wall  away  from  the  separation  line,  then  the 
only  values  of  n not  exceeding  3 for  which  the  vortex  sheet  boundary  conditions 
can  be  satisfied  are  3/2  and  5/2.  The  former  is  associated  with  an  adverse  pres- 
sure gradient  upstream  of  the  separation  line  which  becomes  infinitely  large  as 
the  line  is  approached.  For  n - 5/2  the  pressure  gradient  is  not  dominated  by 
the  local  flow  behaviour.  In  section  5 it  is  shown  that  the  restriction  to 
conical  flow  can  be  removed  without  altering  the  conclusions  and  in  section  6 it 
is  shown  that  the  results  also  hold  for  a vortex  sheet  leaving  a smoothly  curved 
body.  Section  7 includes  a brief  discussion  of  the  results  as  thay  relate  to 
modelling  separated  flow,  both  steady  flow  in  three  dimensions  and  time- 
dependent  plane  flow. 

2 THE  VORTEX  SHEET  LEAVES  THE  BODY  TANGENTIALLY 

In  this  section  it  is  shown  that,  if  vorticity  is  being  shed  from  an 
ordinary  point  of  a separation  line  on  a smooth  body  into  a vortex  sheet  which 
is  embedded  in  a steady,  irrotstional  hcmentrcpic  flow,  then  the  vortex  sheet  is 
tangential  to  the  body  at  that  point. 

By  an  ’ordinary  point’  is  meant  a point  at  which  the  separation  line  has  a 
single,  continuously  varying,  tangent;  and  by  a 'smooth  body'  is  meant  a body 
which,  on  the  separation  line,  has  a single,  continuously  varying,  tangent  plane. 
The  assumption  of  steady  flow  is  primarily  for  simplicity.  'Homentropic'  means 
that  the  entropy  is  the  same  throughout  the  flow  field.  The  assumption  of  homen- 
tropic, irrotational  flow  excludes  the  study  of  vortex  sheets  in  flows  with 
strong  shock  waves,  and  of  vortex  sheets  bounding  closed  bubbles  of  fluid  which 
are  not  penetrated  by  fluid  from  the  free  stream.  These  assumptions  are  likely 
to  be  valid  for  separation  from  a highly-swept  separation  line  on  a stationary 
body  producing  small  disturbances  in  a uniform  stream,  but  it  is  not  necessary 
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to  assume  that  disturbances  are  small  or  that  the  separation  line  is  highly 
swept. 

As  a consequence  of  the  assumption  of  irrotationality  and  the  elimination 
of  closed  bubbles,  the  circulation  of  the  vortex  sheet  can  be  defined.  The 
circulation  along  a path  which  begins  at  a point  on  one  side  of  the  vortex  sheet 
and  ends  at  the  adjacent  point  on  the  opposite  side  of  the  sheet,  lying  otherwise 
in  the  irrotational  fluid,  depends  only  on  the  choice  of  point  on  the  sheet.  If 
the  point  lies  on  the  separation  line  at  the  base  of  the  sheet,  the  circulation 
along  such  a path  is  defined  (with  some  sign  convention)  as  the  circulation  of 
the  sheet  at  that  position  on  the  separation  line.  The  notion  of  the  shedding  of 
vorticity  can  now  be  made  more  precise.  It  is  taken  that  vorticity  is  being  shed 
from  a point  on  the  separation  line  if  and  only  if  uhe  circulation  is  varying 
with  position  on  the  separation  line  at  this  point  and  the  point  is  not  a stagna- 
tion point  of  the  mean  (or  convective)  flow.  An  example  in  which  vorticity  is 

not  being  shed  is  provided  by  the  inviscid  model  of  the  flow  pact  a combination 

of  a lifting  wing  and  a cylindrical  fuselage:  the  plane  vortex  sheet  representing 
the  wake  from  the  wing  lies  along  the  body  side  without  being  fed  from  the  body. 
The  present  results  do  not  relate  to  such  a situation. 

Following  this  explanation,  the  proof  is  very  simple.  We  assume  first  that 
the  sheet  is  not  tangential  to  the  body  and  show  that  vorticity  is  not  being  shed. 

Since  the  body  and  the  separation  line  are  continuously  curved,  there  are 
unique  normals  to  the  body  and  to  the  vortex  sheet  at  the  point  considered. 

These  normals  are  distinct,  since  we  assumed  the  sheet  is  not  tangential  to  the 

body,  and  the  component  of  the  fluid  velocity  along  each  vanishes,  since  both  the 
sheet  and  the  body  are  stream  surfaces  of  a steady  flow.  Hence  the  fluid 
velocity  on  both  sides  of  the  sheet  is  perpendicular  to  both  normals,  ie  it  is 
parallel  to  the  tangent  to  the  separation  line. 

The  second  boundary  condition  on  the  sheet  is  that  the  pressure  is  contin- 
uous across  it.  For  steady,  homentropic,  irrotational  flow,  equality  of  pressure 
implies  equality  of  fluid  speed  (neglecting  body  forces),  so  the  fluid  velocity 
on  the  two  sides  of  the  separation  line  is  either  (a)  zero,  (b)  equal  in  magni- 
tude and  opposite  in  direction,  or  (c)  non-zero  and  identical.  Possibilities 
(a)  and  (b)  both  correspond  to  zero  mean  velocity,  ie  to  a stagnation  point  of 
the  mean  flow,  so  that  vorticity  is  not  being  shed.  This  leaves  possibility  (c) , 
according  to  which  the  fluid  velocity  is  the  same  on  both  sides  of  the  separation 
line*  Hence  the  fluid  particles  which  at  one  instant  lie  on  a path  which 
surrounds  the  vortex  sheet,  beginning  and  ending  at  adjacent  points  on  opposite 
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sides  of  the  separation  line,  at  the  next  instant  also  lie  on  a path  which 
surrounds  the  sheet  and  terminates  at  adjacent  points  a short  distance  downstream 
along  the  separation  line.  However,  the  circulation  along  such  a path  does  not 
change  with  time,  so  the  circulation  of  the  vortex  sheet  is  not  varying  with 
position  along  the  separation  line,  and  so  vorticity  is  not  being  shed. 

This  shows  that  vorticity  is  not  being  shed  if  the  sheet  is  not  tangential 
to  the  body.  On  the  other  hand  we  can  see  that  vorticity  can  be  shed  if  the 
sheet  is  tangential  to  the  body.  It  is  convenient  to  distinguish  between  the  two 
sides  of  the  separation  line,  and  consequently  of  the  vortex  sheet,  referring  to 
the  side  on  which  the  sheet  lies  close  to  the  surface  as  the  downstream  side. 

Then  on  the  downstream  side  the  velocity  vector  must  again  lie  parallel  to  the 
separation  line,  but  on  the  opposite,  upstream,  side  the  only  restriction 
imposed  on  the  direction  of  the  velocity  is  that  it  must  be  parallel  to  the 
conmon  tangent  plane  of  the  sheet  and  the  body.  The  magnitudes  of  the  velocity 
vectors  must  again  be  equal,  but  this  allows  there  to  be  both  a non-zero  com- 
ponent of  the  mean  velocity  normal  to  the  separation  line  and  a non-zero  com- 
ponent of  the  vorticity  parallel  to  the  separation  line.  Vorticity  can  therefore 
be  shed  and  the  circulation  of  the  sheet  can  vary  along  the  separation  line. 

A sketch  of  the  configuration  envisaged  is  shown  in  Fig  1 . The  portions 
of  body  surface  and  vortex  sheet  lying  between  two  parallel  cross-flow  planes 
are  represented,  with  the  vortex  sheet  touching  the  body  along  the  separation 
line.  Streamlines  of  the  inviscid  flow  on  the  wall  and  on  the  sheet  are  shown 
by  solid  lines.  On  the  downstream  side  of  the  separation  line,  these  are  all 
nearly  parallel  to  the  separation  line,  whether  they  lie  on  the  wall  or  on  the 
vortex  sheet.  On  the  upstream  side,  the  inviBcid  streamlines  on  the  wall  cross 
the  separation  line  at  an  angle  and  continue  on  the  vortex  sheet.  The  skin- 
friction  lines,  or  limiting  streamlines  of  the  viscous  flow,  are  shown  on  the 
wall  upstream  of  the  separation  line  as  broken  lines.  These  deviate  inwards 
from  the  curved  inviscid  streamlines,  approach  the  separation  line,  but  do  not 
cross  it.  On  the  downstream  side,  a nearly  parallel  flow  is  envisaged  in  which 
the  skin-friction  lines  on  the  wall  are  essentially  parallel  to  the  inviscid 
streamlines.  The  direction  of  the  skin  friction  thus  varies  continuously  on  the 
wall  and  the  velocity  varies  continuously  on  each  side  of  the  sheet.  The 
boundary  layer  on  the  downstream  side  is  r.ot  obliged  to  separate  in  this  concept 
of  the  flow,  for  which  support  is  provided  in  subsection  4.4. 
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A BASIS  FOR  FURTHER  ANALYSIS 


The  previous  section  showed  that,  under  fairly  general  assumptions,  (steady 
flow,  an  'open*  type  of  separation  and  an  absence  of  strong  shock  waves)  a vortex 
sheet  into  which  vorticity  is  being  shed  from  a smooth  body  is  tangential  to  the 
body  surface  along  the  separation  line.  To  find  out  more  about  the  form  of  the 
vortex  sheet  and  the  flow  field  near  the  separation  line  a more  complete  analytic 
framework  is  needed.  Since,  as  indicated  in  the  Introduction,  flows  with  highly- 
swept  separation  lines  are  of  considerable  interest,  it  is  nstural  to  choose  the 
framework  of  slender-body  theory.  This  means  we  are  concerned  with  a body  of 
slowly  varying  cross-section,  whose  lateral  dimensions  are  small  compared  with 
its  length,  with  its  longitudinal  axis  at  a small  inclination  to  a uniform  stream 
whose  Mach  number  is  not  large  compared  with  unity.  On  such  bodies  extensive 
flow  separations  of  open  type  (is  which  are  penetrated  by  upstream  fluid)  are 
observed,  for  which  a vortex-sheet  model  is  likely  to  be  appropriate.  It  is  to 
such  a model  that  the  results  obtained  will  apply. 

According  to  slender-body  theory,  the  flow  is  described  by  a disturbance 
velocity  potential  which  is  the  sum  of  two  terms.  The  first  is  a function  of  the 
streanwise  coordinate  only,  with  a form  which  depends  on  the  Mach  number  cf  the 
free  stream.  This  function  only  gives  rise  to  a disturbance  velocity  in  the 
streamwise  direction,  and  this  velocity  is  uniform  in  each  cross-flow  plane,  so 
it  emerges  that  it  plays  no  part  in  the  present  analysis.  The  second  term  is  a 
harmonic  function  (a  solution  of  Laplace's  equation)  in  the  cross  flow  plane, 
and  is  therefore  obtainable  as  the  real  part  of  a complex  analytic  function  of 
the  complex  coordinate  in  this  plane.  For  a wholly  attached  flow,  which  is 
completely  irrotational,  this  function  is  determined  (to  within  an  additive 
constant)  in  each  cross-flow  plane  by  the  boundary  conditions  on  the  body  surface 
and  at  infinity.  When  vortex  sheets  are  present  further  conditions  are  required. 
For  instance,  in  a linearized  approach,  where  the  vortex  sheet  lies  in  a fixed 
position  and  represents  the  wake  shed  from  the  trailing  edge  of  a wing,  a 
condition  of  pressure  continuity  is  used  to  determine  its  strength.  For  problems 
of  separated  flow,  where  the  position  as  well  as  the  strength  of  the  vortex 
sheet  is  to  be  found,  the  two  boundary  conditions  that  the  sheet  is  a stream 
surface  and  that  the  pressure  is  continuous  across  it  are  required.  These 
conditions  can  be  written  as  differential  equations  which  express  the  variation 
of  the  position  and  strength  of  the  vortex  sheet  in  the  streamwise  direction  in 
terms  of  the  velocity  field  in  the  cross-flow  plane  (see,  for  example,  the 
Appendix  or  Refs  4 and  5).  For  a sharp-edged  body,  for  which  the  separation  line 
is  known,  these  equations  can  be  integrated  to  trace  the  downstream  development 
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of  a vortex  sheet  specified  at  some  station  upstream^,  It  is  to  be  expected  that 
the  same  will  be  true  for  a smooth  body,  once  the  separation  line  is  specified; 
leaving  the  specification  of  the  separation  line  and  the  establishment  of  an 
initial  vortex  configuration  upstream  as  the  outstanding  problems. 

These  problems  do  not  enter  the  present  analysis,  since  no  attempt  is  made 
to  determine  the  entire  vortex  sheet  configuration.  The  more  limited  aim  of 
further  restricting  the  possible  behaviour  of  a vortex  sheet  near  a separation 
line  can  be  achieved  by  arguments  in  a single  cross-^lov  plane,  like  those  used 
by  Clapworthy  and  Mangier  to  discuss  the  vortex  shi  c leaving  a salient  edge  in 
a conical  flow.  Suppose  that  ve  assume  a shape  for  the  section  of  the  sheet  by 
a cross-flow  plane  and  also  specify  the  variation  with  the  streamwise  coordinate, 
x , of  that  shape  in  the  immediate  neighbourhood  of  the  cross-flow  plane.  Then 
the  component  of  the  velocity,  in  the  cross-flow  plane,  which  is  normal  to  the 
curves  in  which  the  body  and  sheet  intersect  the  plane  in  known  (by 
equation  (A-4) , for  instance)  and  so  is  the  behaviour  of  the  flow  at  infinity 
in  this  plane.  Hence  the  entire  flow  in  this  plane  is  determined.  It  is  helpful 
to  set  out  one  way  of  finding  it  explicitly. 

Let  us  imagine  that  the  entire  region  of  the  cross-flow  plane,  Z - y + ix  , 
exterior  to  the  body  and  the  vortex  sheet  (or  sheets)  is  mapped  conformally  on  to 
the  upper  half  of  a transformed  plane,  Z*  , with  the  point  at  infinity  in  the 
Z*-plane  corresponding  to  the  point  at  infinity  in  the  Z-plane  by  an  analytic 
function 

Z " F(Z*)  . 

In  this  mapping  the  surface  of  the  body  and  both  surfaces  of  the  sheet(s)  become 
part  of  the  real  axis  in  the  Z*-plane.  We  now  construct  the  flow  field  in  this 
transformed  plane.  A source  distribution  is  introduced  along  the  real  axis 
whose  strength  per  unit  length  is  proportional  to  the  product  of  the  normal 
velocity  in  the  cross-flow  plane  and  the  modulus  of  the  mapping: 

q(y*)  - 2v*(y*)  - 2vn(Z) |F'(y*)|  , (1) 

where  Z is  the  point  on  the  configuration  corresponding  to  the  point  Z*  « y*  , 
real,  ie  Z »■  F(y*)  . Then  the  complex  potential  is  given  by 

ao 

. VZ*  + f q(y*>  In  (Z*  - y*)dy*  , 


W 


(2) 


1C 


where  V is  the  speed  of  the  flow  at  infinity  in  the  Z*-plane.  The  component  of 
velocity,  in  the  cross-flow  plane,  which  is  tangential  to  the  curves  in  which  the 
body  and  sheet  intersect  the  plane  is  then  given  by 


vt<z0)  - v£(y$)/l*' (y$> I » (3) 

where  ZQ  - F(y*),  y*  real,  and  where  v*  , the  corresponding  velocity  component 
in  the  transformed  plane,  is  given  by  the  derivative  of  (2): 


v*(y§) 


V 


+ _L  [ q(y*>dy* 
2n  J y§  - y* 
—00 


(4) 


The  integral  in  (4)  is  to  be  interpreted  as  a Cauchy  principal  value.  From  the 
values  of  v at  points  on  the  sheet,  the  right-hand  side  of  equation  (A-10), 


3A» 

dX 


. . 3z  . , ml 

COS  ill  + — Sin  ill 77 

y 3x  r U 


(A-10) 


expressing  the  condition  of  pressure  continuity,  can  be  found.  So,  too,  can 
A*  . If  this  determination  can  be  carried  through  in  terms  of  parameters  whose 
variation  with  the  streamwise  coordinate,  x , is  known,  the  differentiation 
on  the  left  of  (A-10)  can  be  performed  and  the  consistency  of  the  assumed  shape 
of  the  sheet  with  the  pressure  condition  can  be  checked. 

This  procedure  is  not  a practical  way  to  investigate  the  global  solution, 
but  it  can  be  applied  to  obtain  local  information  at  points  of  the  flow  field 
at  which  the  velocity  is  singular.  To  make  this  idea  more  precise^  we  recall 
that  a complex  function  of  a complex  variable  is  regular  at  a point  if,  and  only 
if,  it  can  be  expanded  in  a power  series  about  that  point,  the  power  series  hav- 
ing a non-zero  radius  of  convergence.  If  a function  is  not  regular  at  a point, 
it  is  singular  there.  The  same  definition  can  be  applied  to  real  functions  of  a 
real  variable,  with  the  power  series  being  convergent  in  a non-vanishing  interval 
centred  on  the  point.  (Note  that  a function  which  is  singular  in  this  sense  may 
be  bounded,  like  the  Heaviside  step  function,  and  even  differentiable,  like 
exp  (-  1/x  ) , both  at  x ■ 0 .)  There  is  clearly  a sense  in  which  a singular 
function  is  unchanged  when  a regular  function  is  added  to  it;  we  shall  say  that 
its  singular  behaviour  is  unchanged.  On  a more  formal  level  we  could  define 
equivalence  classes  of  functions,  assigning  functions  to  the  same  class  if  they 
differ  by  a regular  function;  but  this  refinement  is  not  necessary.  Note  that 
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the  singular  behaviour  of  In  sin  x and  In  x at  x ■ 0 is  the  same,  but  that 

2 

of  2 In  x is  different;  nor  is  the  singular  behaviour  of  (In  sin  x)  the  same 
2 

as  that  of  (In  x) 

We  can  now  prove  a theorem  which  is  fundamental  to  the  analysis  which 
follows,  viz;  If  v*  is  given  by  (4),  then  its  singular  behaviour  at  a point 
P is  uniquely  determined  by  the  singular  behaviour  of  q (or  v*)  at  P . 

The  proof  follows. 

We  choose  the  origin  at  P and  simplify  the  notation,  writing  (4)  as 

00 

v(o  - v ♦ JL  J ilaiii  . (4') 


We  wish  to  prove  that,  if  q j (n)  and  q2(n)  have  the  same  singular  behaviour 
at  n ■ 0 and  Vj(£)  and  V2(C)  are  related  to  them  by  (4'),  then  Vj  and 
V2  have  the  same  singular  behaviour  at  £ » 0 . Let  q(ri)  * qj(n)  ~ qj^)  and 
v(£)  * Vj(5)  - V2(S)  • Since  (4*)  is  linear,  q and  v are  also  related  by  it. 
Then,  by  the  definition  above  of  singular  behaviour,  we  must  prove  that  the 
regularity  of  q at  n * 0 implies  the  regularity  of  v at  £ - 0 . If  q 
is  regular  at  the  origin  we  can  write,  for  some  e > 0 , 


q(n) 


n»0 
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for  |n | < e . 

Now  suppose  j?)  < e and  rewrite  (4*)  as 

*"«>  ■ **♦  Tt“**  / E v. . 

-®  e n"0 

"here  * / H’V  ■ 

-e 

Writing  n”  ■ Cn”  1 “ (£  “ n)nn  1 , we  see  that 

I - Cl  . - 2en/n  for  n > 1 , 

n n 1 


while 
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Hence  Iq  ■ £n  In  + a polynomial  in  5 , and  so  I 

power  series  in  K for  ] £ j < e 


In  the  integrals  over 
the  expression  for  2irv(5)  , we  can  expand  (5  - n)  * as 
|e|  < e < n , giving 


can  be  expanded  as  a 
infinite  intervals  in 
a power  series,  for 


f q(n)dn 

J 5 - n 


-1*1 


g(n)dn 

n+1 

n 


and  a similar  expression.  Hence  v(£)  is  regular  at  £ ■ 0 and  the  proof  is 
complete. 


The  significance  of  this  theorem  is  that  singular  behaviour  is  likely  to 

arise  on  the  two  sides  of  the  separation  line.  This  may  coca  from  the  shape  of 

the  cross-section  of  the  body  and  sheet,  through  the  mapping  function  F in 

(1)  and  (3),  or  from  the  streamwise  rate  of  change  of  the  cross-section,  through 

the  normal  velocity  v in  (1).  Moreover,  the  behaviour  of  q or  v*  is 

n ti 

locally  determined,  through  (1).  Hence  the  theorem  tells  us  that  the  singular 
behaviour  of  v*  on  the  two  sides  of  the  separation  line  is  locally  determined. 
This  can  be  used  to  provide  information  on  the  local  behaviour  of  v , through 
(3),  and  the  consistency  of  this  behaviour  with  equation  (A-1Q)  can  be  checked. 
This  provides  the  main  result,  that  a number  of  plausible  forms  for  the  vortex 
sheet  must  be  rejected. 


The  last  preliminary  question  to  be  discussed  before  embarking  on  the 
analysis  concerns  the  mapping  function,  F(Z*)  . The  argument  above  only 
involves  properties  of  F in  the  neighbourhoods  of  the  two  sides  of  the  separa- 
tion line.  It  is  therefore  sufficient  to  construct  a function  which  maps  the 
region  outside  the  body  and  a truncated  vortex  sheet  on  to  the  upper  half  of  the 
Z*  plane,  provided  that,  near  the  separation  line,  the  truncated  sheet  has  the 
local  behaviour  which  is  to  be  investigated. 

The  next  section  begins  with  the  construction  of  an  appropriate  mapping. 

The  argument  is  then  carried  through  for  a conical  sheet  on  a plane  wall,  a case 
which  brings  out  the  essential  features  without  unnecessary  complications.  In 
sections  5 and  6 the  limitations  to  conical  flow  and  plane  walls  are  removed. 

The  reader  who  is  prepared  to  take  the  details  on  trust  may  proceed  to  section  7. 

4 THE  CONICAL  SHEET  ON  A PLANE  WALL 

4 . 1 The  sheet  shape  and  mapping  function 

In  this  section  and  the  next,  the  body  shape  considered  is  a plane  wall 
which  extends  laterally  to  infinity  and  is  parallel  to  the  undisturbed  flow.  It 
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does  not,  of  itself,  disturb  this  flow,  so  slender  body  theory  is  applicable, 
provided  the  vortex  sheet  is  slender.  Only  highly-swept  separation  lines  are 
considered  and  it  is  assumed  that  a vortex  sheet  springing  from  a highly-swept 
separation  line  is  slender. 

The  shape  of  the  sheet  is  described  in  a cross-flow  plane  which  is  normal 
to  the  separation  line.  In  this  section,  the  sheet  is  assumed  to  be  conical, 
ie  to  be  generated  by  straight  lines  through  an  apex,  so  that  its  entire  shape 
is  defined  by  its  section  by  the  cross-flow  plane,  together  with  the  distance 
of  that  plane  from  the  apex.  In  this  case  it  is  convenient  to  choose  the  origin 
0 , of  the  right-handed  rectangular  Cartesian  system  of  coordinates  at  the  apex, 
with  Ox  downstream  along  the  separation  line,  0y2  to  starboard  and  Oz^  normal  to 
the  plane  wall  and  into  the  fluid,  as  sketched  in  Fig  2.  We  suppose  y2  < 0 is 
the  upstream  side  of  the  separation  line  on  the  wall,  *2  ■ 0 . Then,  near  the 
separation  line,  ie  for  small  values  of  y2  and  z2  , the  section  of  the  sheet 
by  the  cross-flow  plane  x ■ const  is  given  by 


*2  ~ uy2  * 


n > 1 . 


(5) 


In  this  Report  the  tilde  is  used  to  indicate  asymptotic  equality,  ie 


a(t)  ~ b(t)  i'.f  and  only  if  lim 


ft. 


a(t) 

b(t) 


The  appropriate  variable,  t , and  limiting  value,  t^  , is  usually  obvious  from 
the  context.  When  an  indication  of  difference  between  the  left  and  right  hand 
sides  is  desired,  the  0 notation  is  used: 


a(t)  - b(t)  + 0(f(t)) 


if  and  only  if  there  exist  A and  e such  that 

| a(t)  - b(t) | « A j f (t) | for  |t  - tQ|  < e . 

As  explained  in  the  previous  section,  it  is  sufficient  to  replace  the  sheet  by 
any  arc  which  has  the  same  behaviour  near  the  separation  line.  We  obtain  an  arc 
in  the  plane  of  Z2  - y2  + iz2  with  the  behaviour  (5)  by  taking  a circular  arc 
in  a plane  Z^  , related  to  Z^  by 
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m > 0 , (6) 

where  a is  the  diameter  of  the  circular  arc.  To  obtain  the  required  plane 
wall  in  the  plane,  we  introduce  a wall  with  a corner  at  the  origin  in  the 
plane,  as  shown  in  Fig  3.  In  the  Z^  plane,  the  fluid  occupies  the  sector 

0 < e3  < v/m  , 

where  Z^  ■ r3e*®3  • 

m must  be  at  least  $ if  the  Z^  plane  is  not  to  be  covered  twice,  a complica- 
tion we  choose  to  avoid.  The  half-lines  0^  ■ 0 and  0^  • ir/m  represent  the 
wall,  and  the  vortex  sheet  is  represented  by  a circular  arc  of  radius  a/2  , 
touching  the  wall  0^  - 0 at  the  origin.  A general  point,  P , on  the  sheet  is 
given  by 

r^  “ a sin  03  , 0 < 0j  < Bq  (7) 


where  the  suffix  C denotes  the  value  of  »3  at  the  end  of  the  arc.  Suffixes 
‘U  and  d denote  the  upstream  and  downstream  sides  of  the  sheet  at  P . A is 
the  point  at  infinity,  B and  D represent  the  upstream  and  downstream  sides  of 
the  separation  line. 


Corresponding  points  are  denoted  by  the  same  letters  in  the  plane. 

Under  (6)  the  sector  becomes  the  upper  half -plane  0 4 0O  < it  , where 

i02  ^ 

Z2  " r2e  ’ *-8  t*ie  *eal  axis  and  the  equation  of  the  vortex  sheet  is, 

by  (6)  and  (7) 


. m 

a sin  9, 


a sin 


_«!i 
m 


(8) 


For  small  values  of  0^  , (8)  yields 


_2_ 

a 

f2 

a 


COS  0O 

a 1 


■ — — sin  9„  ~ 


,m+l 


m 

m 


( 1 +m)  /m 


(9) 


The  required  beheviour  (5),  is  therefore  obtained  with 


m ■ 


and 


Thus  the  order  of  contact  increases  as  m falls.  The  curvature  of  the  sheet  at 
the  separation  line  is  zero  for  m < 1,  n > 2 ; finite  for  m ■ 1,  n ■ 2 ; and 
infinite  for  m > 1,  1 < n < 2 . 


The  region  occupied  by  the  fluid  in  the  plane  is  bounded  by 
lines  and  circular  arcs,  all  with  a common  point.  It  can  therefore  be 
on  to  a region  bounded  by  straight  lines  by  an  inversion  which  carries 
point  to  infinity: 


straight 
napped 
the  common 

(10) 


The  configuration  in  the  plane  is  shown  in  Fig  4.  The  half-lines  BA  and  DA 
remain  half-lines,  though  the  origin  is  interchanged  with  the  point  at  infinity. 
The  circular  arc  becomes  a third  half-line,  parallel  to  the  reel  axis,  since  the 
arc  touches  the  real  axis  at  the  origin.  If  Z^  lies  on  the  sheet, 

Zj  - a sin  e^e*^  , by  (7),  and 


Z^  ■ a(cot  63  - i)  » 


by  (10).  A convenient  relation  is  obtained  from  (6)  and  (10): 


a 


(ID 


Since  the  region  occupied  by  the  fluid  in  the  Z^  plane  is  bounded  by 
straight  lines,  it  can  be  mapped  on  to  the  upper  half  of  the  Z^  plane  by  a 
Schwarz-Christoffel  transformation.  If  A maps  into  the  point  at  infinity  and 
B,  C and  D into  - b,  c and  d , where  b < c < d , as  in  Fig  4,  we  find  the 
mapping  function  is  defined  by 


1 Z5  ~ c / a 
m d - Z5  ^Z5  - bj 


(12) 


where  a disposable  constant  has  been  chosen  so  that  Z2  ~ Z^  at  infinity.  The 
rules  for  constructing  (12)  are  given  in  standard  texts  which  deal  with  conformal 
mapping. 
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It  is  not  easy  to  integrate  (12)  for  general  values  of  m . Fortunately, 
it  is  sufficient  for  our  purpose  to  consider  the  behaviour  of  the  transformation 
near  B and  near  D.  Let  us  first  consider  the  neighbourhood  of  B and  vrite: 


Z5  - b ♦ ac  , 


C small  . 


(13) 


Introducing  this  into  (12),  we  have 


<*Z4 

dC 


mC 


(m+O/a 


(1  - k2?  - k3C  ♦ ...)  , 


(14) 


where  the  omitted  terms  are  regular  and  0(cJ>  and 


|t  m jLJL.il  a 

K1  d - b » *2 


a(d  •*  c) 

(c  - V>(d  -b7 


and 


k3  " 


az(d  - c) 


k2  * 


(c  - b)(d  - b)‘ 
(14a) 


In  integrating  (14),  it  is  necessary  to  consider  separately  the  cases 


m » J , i < m < 1 , m - 1 and  m > 1 . 


(a)  m - i . Equation  (14>  becomes 


therefore 


Z4  - ak 


3T  • 2*ki(-7+7*T*0<0) 

*(?" 


+ 2kj  In  C + const  + 0(c) 


-) 


Hence,  by  equation  (U)» 


a 


-Jj-  (l  + k2c  " k3c2  ln  5 ♦ AC2  ♦ 0(c3  In  C>) 
kl 


05) 


where  A is  a complex  constant,  so  far  arbitrary.  Now,  on  BA,  C i«  real  and 
negative,  so  ln  c ■ In  |c|  ♦ iw  , while  Z.  must  be  real.  Hence  the  imaginary 
psirt  of  A must  cancel  the  imaginary  part  of  the  preceding  term  in  (15): 
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A - k3(Ar  + iw)  , 


where  Af  is  an  arbitrary  real  conatant.  Also,  on  BC,  c ia  real 
ao,  by  (IS)  and  (16) 

• r *»  > 

1,3 


and 


1 

Comparing  theae  expressions  with  (S)  and  (9),  we  see  we  need 

2irk]k3  « 1 or  2ira^(d  - c)  » (d  - b)^  . 


(b)  | < m < 1 . Equation  (14)  becomes 


therefore 


therefore 


““4  “1  / 1 k2  _(nr-l)/mA 

d r “r7WF«;  + 7i7s+  0(c  >) 


z,  - 


i 

?l/m 


('  - r~  ♦ w17" ♦ o<c2)j 


Z0  _ / mk„C 

-±  . JL  h + — - 

. k»  \ 1 -■ 


l*c1/n  ♦ 0(c2) 


Again,  on  BA  we  require  Z„  real,  with  C ■ |c|elir  . Hence 


A - |A|e‘i”/“ 


Also,  for  C real  and  positive,  corresponding  to  BC,  we  require 


and  so 


■*iz4}  - - « » 

- a - akj/{A}  . 


06) 

and  positive. 


(17) 


(18) 

(19) 


Hence,  by  (19), 
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and,  by  (18) 


t ■ ^ ('  * r=r ' “(cot « ' ¥'%  * 0('2))  • 


(c)  m - 1 . Equation  (14)  becomes 


dZ, 


dC 


- ak 


(- Vt*0(1>) 


therafora 


~~  ^1  + *2c  In  C ” AC  + 0(C2)^  » 


where  A is  an  arbitrary  complex  constant.  Hence,  by  (11), 


T - 


k2C  In  C + AC  ♦ 0(c  In 


»)  . 


As  in  (a)  above,  c real  and  negative  corresponds  to  Z2  real  and  so 


A - k2(Ar  ♦ is)  , 


where  Ar  is  real . 


Also,  £or  C real  and  positive. 


*k. 


2 r2 


‘I 


end  so,  by  comparison  with  (9): 


irkjk2  - 1 


or 


sa(d  - c)  " (d  - bj  . 


(d)  m > 1 . The  expression  for  Z . in  (b)  above  becomes 
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therefore 


The  same  arguments 


h . JL  (,  - „*•/»  + 0(5,  . 

a . m \ t 

K1 

as  in  (b)  lead  to  the  same  value  of  A , 


so  that 


a 


(24) 


We  have  now,  in  equations  (13)  to  (24),  expressions  describing  the  behaviour  of 
the  mapping  near  B , for  all  relevant  values  of  m . 


Near  D the  form  of  the  expression  does  not  depend  on  m . Let 


Z 


5 


d + au  , 


u small  . (25) 


Then,  by  equation  (12) 


where  a 
u 

Therefore 

Now,  for  a real  and  positive,  corresponding  to  DA,  is  real,  so  A is 

real.  Let  B ■ eA  , so  now  B is  an  arbitrary  positive  constant.  For  io  real 
and  negative,  corresponding  to  DC,  /{z^f  ■ -a  , and  so 


- a 


(m+I)/m  d _ 


m 


(27) 


- - ¥ ln  B"  (j + °fcrs:)) 


Hence  (26)  becomes 
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and,  by  (li) 


t ’ (r-dhc)"  ('  ♦ °(i she))  • <“> 

4.2  T.ocal  behaviour  of  velocity  field 

The  first  step  along  the  line  of  approach  described  in  section  3 is  to 
obtain  the  normal  velocity  on  the  sheet  in  the  cross-flow  plane,  Z ^ . In  the 
Appendix,  the  boundary  condition  that  the  sheet  is  a stream  surface  ia  used  to 
provide  expressions  for  this  normal  velocity  component.  There,  as  here,  the 
cross-flow  plane  is  normal  to  the  axis  of  x , which  is  inclined  at  a small  angle 
to  the  undisturbed  stream.  If  the  sheet  is  conical,  vith  apex  at  the  origin, 
equation  (A-5)  becomes  (see  Fig  3 for  notation) 


v 


where  the  sense  of  the  normal  is  given  by  a positive  rotation  of  w/2  from  the 
direction  of  6 increasing  along  the  curve,  ie  into  the  fluid  on  BC  and  into 
the  sheet  on  DC.  Note  that  r,  8,  $ and  v^  all  refer  to  the  Z^-plane,  but 
suffixes  are  omitted  for  simplicity . Equation  (8)  relates  r to  6 : 

t - . . 

If  o is  the  arc  length  along  the  sheet  in  the  Z2~plane,  measured  from  the 
separation  line. 


» 


and  sin  4 ■ rdd/do  , so  the  boundary  condition  becomes 


- 1 .in"*1-'6 


(S) 


si*r>M 


(29) 


To  proceed  to  the  next  step,  finding  the  normal  velocity  v*  (or  the 
source  strength  q)  in  the  Z^-plane,  through  equation  (1),  we  need  the  modulus 
of  the  mapping  at  points  on  the  sheet.  By  (11)  and  (12) 


C - z, 
d - Z. 


t&l 


(m+l)/m 


(30) 


We  need  the  local  behaviour  of  tha  modulus  of  this  on  tha  real  axis  near  P and  D, 
la  for  C and  u>  small  and  real.  Since,  on  the  sheet,  do/d - IdZj/dZ^I  , 
it  is  convenient  to  relate  a to  C and  u at  the  same  time.  As  in  the 
previous  sub-section,  we  deal  first  with  the  neighbourhood  of  B,  for  the  differ- 
ent values  of  m . By  (13)  and  (30) 


. . (mH)/m 

«-•>-« (M 

d - b - ac  \mt) 


(31) 


(a)  m ■ | . Introducing  (15)  into  (31),  and  taking  the  modulus  for  C real, 
ve  find 


-J  (>  ♦ 2k2C  - 3k3c2  In  | C | ♦ 0(C2)) 
kl 


Integrating  (32)  for  C > 0 gives 


I " "V  (*  * k2C  " k3c2  ln  C * 0(c2>)  ' 

kl 


(32) 


(33) 


(b)  1 < m < 1 . Introducing  (20)  into  (31)  and  taking  the  modulus  for  c real 

and  positive,  we  have: 


2 

•a. 


-["T1 

•T\ 


c - 


wu  I u 

in  1 


Integrating  (34)  for  C > 0 gives 


£ 

a 


(34) 


(35) 


(c)  m ■ 1 . Introducing  (21)  into  (31)  and  taking  the  modulus  for  C real, 
ve  find: 
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dZ, 


dZ, 


£j-  (»  - 2k2C  In  |c|  + 0(C))  . 


Integrating  (36)  for  C > 0 gives 


5“  - ^ (j  “ ? j5  ln  5 + 0(C))  * 


(d)  m > 1 . Introducing  (24)  into  (31)  and  taking  the  modulus  for  c 
and  positive,  we  have: 


dZ. 


dZ. 


. -L  (\  - aJ iJ. 

*?V  k' 


cot  J |c1/m|  ♦ 0(C 


Integrating  (38)  for  c > 0 gives 


S ■ i (•-£-*  i i«,mi  * •«.***>)  • 

*1  1 

The  form  of  the  leading  tanas  in  (32)  to  (39)  is  independent  of  m , but 
higher  order  terms  differ  significantly. 

for  the  neighbourhood  of  0,  introducing  (25)  into  (30)  gives 


dZ2 

dz! 


c - d - aw 


- aw 


L h \ 

Id  - b - aw J 


(m+l)/m 


By  (27)  and  (28)  this  reduces  to 


dZ2 

dZl 


• £ {- ttisT' ° * 0M)  ■ 


On  the  sheet  near  D,  w is  real  and  negative,  so  that 


ln  Bw  ■ ln  B|w|  + iff  ■ - ff(l  - ifi)/5  , 


where 


6 


ln  B|w| 

is  a convenient  small  positive  quantity.  The  modulus  of  (40)  becomes 


dZ. 


dZ, 


m6 


m+1 


L|<1  ♦ j2) <■**'>« 


(1  + 0(w))  for  w < 0 


and 


(36) 

(37) 

real 

(38) 

(39) 

the 


(40) 

(41) 

(42) 

(43) 
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(1  + 0(w» 


for  w > 0 . (44) 


Not*  that , by  (42),  u ia  exponentially  small  in  relation  to  4 . Also  by  (42), 
2 

dfi/d u - 6 /iru  and  so,  for  u < 0 : 


"2 


|dZ,5|  |do>  |Jd« 


(I  + O 


(1  ♦ 0(w)> 


o - a«“(l  + 0(62))  . 


We  can  now  find  the  component  v*  of  the  velocity  normal  to  the  real  axis 

n 

in  the  Z^-plane,  close  to  B and  D.  By  (t)  and  (29)  we  have,  on  BC  and  DC, 


-£  . - 5.  (1) 

U x U ) 


(m+ 1 ) /ml  dZ 


~ (I  + 0 (o2/®))  , 


with  the  normal  into  the  fluid  on  BC  and  out  of  it  on  DC.  On  BA  and  DA,  which 
correspond  to  the  body,  (A~5)  and  (I)  show  that  v*  ■ r » 0 . Introducing  into 

XX  XX 

(46)  the  leading  terms  from  (32)  to  (39),  (43)  and  (45),  we  find  for  the 
velocity  component  w*  normal  to  the  real  axis  end  into  the  fluid  in  the 


Z*-plane: 

on  BA 

w*  - 

o 

■ 

A 

| 

on  BC 

w*  - 

V*  - - SE 1 

n x k 

$ 

on  DC 

V*  » 

a.  maU 

- v*  * ■ ■ tt  1 
n irx|o)| 

| 

on  DA 

w*  - 

v*  «•  0 . 

(m+l)/m 


(l  + 0(c1/o,c  In  o) 


(1  + 0(5  j) 


We  now  wish  to  find  the  tangential  velocity,  v*  , on  the  real  axis  in  the 
Z^-plan*  near  B and  D.  The  combination  v*  - iw*  is  an  analytic  function  of 
Zj  . If  (v*  - iw*)g  has  the  singular  behaviour  of  v*  - iw*  at  B or  D,  then 


we  can  write 


v*  - iw*  » (v*  - xw*)s  + (v*  - iw*)^  , 
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where  the  second  term  is  regular.  This  term  is  not  determined  by  local  arguments 
end  we  writ**  it  as 

00 

(v*  - iw*)R  - £ Ab 

n-0  n^O 

The  conditions  (47)  on  BA  and  DA.  imply  that  Ag  and  Bq  are  real.  The  singular 
behaviour  cf  v - iw  at  B and  D is  determined  by  . .e.  singular  behaviour  of  w*  , 
by  the  theorem  proved  in  section  3.  We  therefore  proceed  by  constructing  analy- 
tic functions  whose  imaginary  parts  behave  like  w*  near  B and  D.  As  before, 

we  start  with  B for  the  familiar  four  eases. 

(a)  m - | . By  (47),  for  c real, 


(near  B) 


or 


I 


B*. 

w 

XI 


(near  D) 


(49) 


w*  ~ ^ 


x k2 
1 


for  C > 0 


for  C > 0 


3 

Consider  the  function  f(c)  " C (iff  - In  5)  . For  C real  and  ne  ative, 


and  for  C real  and  positive, 


" 0 , 

» ire3  • 


Hence 


(v*  - iw)g  - (iff  - In  c) 

kl 


(50) 


(b)  J < m < 1 . By  (47) , for  c real 


w*  ~ r 


- ®!i  4 

x . 2m+l 


V.  k 


for  C < 0 


for  c > 0 , 


where  2 < y ■ (m  + l)/m  < 3 . For  C ■ pe  , define  C ■ P ® • ***** 

consider  the  function 


f(C)  ■ « ^'^c'Vsin  yir  ■ CW(cot  yir  - i)  . 


For  ( reel  and  negative,  s?  ■ tt  and  J*  jf } **  0 . For  £ real  and  positive, 
d - 0 and  J{t\  « - Cw  , real.  Hence 


(v*  - iw*)g  ~ “ (i  " cot  -) 


since  cot  yir  - cot  — • 

n 


(c)  n - 1 . By  (47),  for  C real 


v*  ~ C 


_ aU 

” x ? 
1 


Proceeding  as  in  (a)  above,  we  find 


(v*  - iw*)_  ~ *sr  (i*  -lat)  ■ 


(51 


for  C < 0 


for  C > 0 


(5 


(d)  m > 1 . The  argument  in  (b)  applies,  with  1 < y < 2 , leading  again  to 
(51). 


Now  in  the  neighbourhood  of  D,  (47)  gives,  for  u real, 


- 252.  * o«2» 

irxw 


w*  • < 


for  to  < 0 


for  (a  > 0 


Consider  the  function  f^Oa)  ■ ("  w/ln  Bu)  /w  . For  w < 0 , 


In  Bu  ■ — tt ( 1 - ifi)/6 


by  (41),  and  so 
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fj^w)  - 6k(l  - i<5)~k/w 

and 

./{fjJ  - kfik+1  (1  + 0(62))/u>  . 

For  u > 0,  ^{f^}  * 0 • Hence 

^ - i«*>s  ■ rsrry=  (-  TZTsf*"  (' + °(-  TT^i)  • «*» 

Equations  (46)  and  (49),  with  the  appropriate  equation  from  (50)  to  (53)  give 
the  complex  conjugate  velocity  in  the  Z,. -plane  near  B and  D. 

We  can  now  complete  the  calculation  by  using  these  expressions  to  find  the 
tangential  velocity,  v , on  the  sheet  in  the  Z2*-plane,  using  a form  of 
equation  (3) : 


v 


t 


vt 


(54) 


in  which  v*  ■ v*  on  BC  and  v*  « -v*  on  DC.  Again  we  treat  first  the  four 
cases  which  arise  on  BC,  where  $ is  real  and  positive. 

(a)  m - J . By  (32) , (48) , (49)  and  (50) 

Aq  + AjC  + A252  - -2iL  53  In  ; + ... 

TTXkj 

V m i ■ ■ ■ ■■■■  — ■■  i i i—  ■ 1 1 ■ »!■ — ■ 

C -j  (l  + 2k2C  - 3k3C2  In  C + 0(c2)) 
kl 

- kj(A0  + A|5  + 3A0k3C2  In  C + 0(i;2))  , 

where  Aj  is  another  real  constant.  Inverting  (33)  to  give 

? - kj  f (1  + 0(a))  , 

1 a 

and  assuming  for  the  moment  that  Aq  ^ 0 , we  can  write,  using  (17), 

vt  - k\kjt  * Aj'o  + £ (|)  In  a * Ofo2)}  , (55) 


where  A" 


■ruuKittf'i&iti 


is  another  undetermined  real  constant 
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(b)  i < m < 1 . By  (34),  (48),  (49)  and  (51): 


Aq  ♦ AjC  + A2c  - 


aU 


xk 


2m+l 


cot  - t<"*l>/B  ♦ .. 


^ ^ * <*«*>) 


‘i  * A! 


t . J-f=  cot  I t1'”  * 0(C2) 


') 


Using  the  inverse  of  (35)  and  assuming  Aq  =£  0 , we  have 

vt  " k?o(j  + A'i°  + 0 + “>  cot  s (£)  + 0(°2>)  * 

where  Aj  and  A"  are  again  undetermined  real  constants. 

(c)  m - 1 . By  (36),  (48),  (49)  and  ^52) 


(56) 


aU  2 

Ay.  + A.C ^ C In  C ♦ ... 

U 1 irxkj 

Vfc  k7  0 ~ 2V  In  C + 0(O) 

a 

Using  the  inverse  of  (37)  with  (23),  and  still  assuming  Aq  # 0 , we  have: 

vt  - kiAo('  * I f ln  ° * 0(<’>)  • <57) 

(d)  m > 1 . Using  (38),  (48),  (49)  and  (51)  we  obtain  a result  like  (56): 

vt  - *®A0^  + (1  ♦ m)  cot  ^ + 0(a)^  . (58) 

It  is  significant  that  the  dominant  terms  in  v arise  in  each  case  from  the 
regular  behaviour  of  v , for  which  the  real  axis  in  the  Z^-plane  is  a stream- 
line, ie  the  dominant  terms  in  v near  B are  of  the  sane  form  as  they  would 
be  if  the  sheet  were  a streamline  of  the  cross-flow. 

We  now  consider  the  behaviour  of  v on  DC,  near  D,  where  a is  real  and 
negative.  By  (53)  and  (41) 
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VS 


maU6 
(2m  + 


2nr*-l 
1 Virno 


0 


0(62)) 


By  (48)  and  (49),  v*  does  not  contribute  to  v*  to  this  order,  and  so,  by 
(43)  and  (54): 


maUfi2m+1 
(2m  + l)irxw 


w 


m6 


m+1 


0 + 0(6  ;) 


We  can  write  this,  by  (45),  as 


v 


t 


Uq 

(2m  + 1 )x 


O + o (o2/my, 


(59) 


In  contrast  to  the  situation  near  B,  the  dominant  terms  in  v£  near  D arise 
from  the  singular  behaviour  of  v , so  that  they  are  uniquely  determined  by  the 
shape  of  the  vortex  sheet  near  D. 


Having  determined  vt  near  D,  we  can  justify  the  assumption  that  Aq  is 
non-zero.  Equation  (59)  shows  that  vfc  vanishes  at  D,  while  (55),  (56),  (57) 
and  (58)  show  that  vt  also  vanishes  at  B if  Aq  « 0 . In  terms  of  the  dis- 
cussion in  section  2,  this  would  mean  that  no  vorticity  is  being  Bhed  from  the 
separation  line.  In  general  this  would  mean  that  the  point  on  the  separation 
line  was  exceptional,  but  for  the  conical  flow  postulated  in  this  section  it 
would  mean  that  no  vorticity  is  being  shed  anywhere  on  the  separation  line. 
Hence  AQ  # 0 . 

The  convective  velocity  + vtl  at  the  separation  line  has  the 

same  sign  as  An  , so  we  conclude  that  A_  > 0 if  vorticity  is  being  shed. 

U v/ 

4.3  Consistency  of  velocity  field  with  continuity  of  pressure 


For  a conical  field,  the  condition  that  the  pressure  is  continuous  across 
the  vortex  sheet,  (A- 10),  becomes 


A* 


(60) 


since  A#,  y and  z are  proportional  to  x , y * r cos  0,  z - r sin  0 and 
♦ • ♦ - 0 • We  wish  to  consider  whether  this  relation  is  consistent  with  the 
tangential  velocity  components  found  in  subsection  4.2  in  the  neighbourhood  of 
the  separation  line.  By  steps  like  those  leading  to  (29),  we  find 


r cos  $ 


o(l  + 0 <0*'“)) 


(61) 


• ui  v u 

- a sin  — cos  — - 


Examining  equations  (55)  to  (58),  we  see  that  the  tangential  velocity  on  the 
sheet  near  B can  be  written  as 


vtlBC  * *0^  * A"a  * <m  * 0(fJ,/m£(o)  + o(o1/mf(a)j 


(62) 


where  A^  ■ A^(c  - b)m/(d  - b)“  # 0 and  A"  are  undetermined  real  constants; 

j—  In  o for  a • | or  ! 

f (a)  - < 

icot  1 
l m 


ko 


(63) 


for  m > | and  ^ 1 ; 


and  the  notation  o(g(o))  stands  for  terms  of  higher  order  than  g(o)  . On  the 
sheet  near  D,  the  tangential  velocity  is  given  by  (59)  as 


1 1 DC  ■ 73.  °°t>»  (l  * • 


2/m, 


(64) 


All  the  terms  in  (60)  can  now  be  represented  for  small  values  of  a . By 
(61),  (62)  and  (63): 


r cos  ♦ - jjvt  - r 


COg  * ~ 2U  (VJbC  + vJdc) 

- - ^0  + Cjo  + (m  + l)(^y ^f (a)  + o(ol/af(o))J  , 


(65) 


where  Cj  is  another  arbitrary  real  constant.  Similarly 


Av^ 


Vt ' BC  “ VOC 


- AJ^  + C2o  + (m  ♦ ,)(f//nf(o)  + o(o,/mf(o))j  , 

where  C2  is  a constant  related  to  C}  . Integrating  (66)  gives 
a 

/ 4v 


(66) 


A4>  - r ■ 


.do 


" An(J  + ic?°2  + 


(uH-l)/m 


f (o)  + o (o 


(mfl )/m 


f (a)) 


) ■ 


(67) 
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«here  r - ~ *D  tot*^  circulation  in  ^he  sheet. 

When  (63),  (66)  end  (67)  are  introduced  into  (60),  it  becomes 


T ♦ 


A- for  * |C2c2  . »(£) 

-.2  , 


(m+l)/m 


f(o)  ♦ o(o 


(° 


,(m+ 1 ) /m 


*<">)) 


■ ” ^ ♦ (Cj  ♦ 0^)0  + 2(m  + 1)^~j  f(o)  + o(o^m£(o)^ 


There  is  no  term  on  the  left-hand  side  of  (68)  to  balance  the  term 

fm 


,J  ,o\U\.  x 

— u — fe;  £(a) 


(68) 


on  the  right,  so  this  term  must  vanish.  Hence,  since  A£  # 0,  f(o)  ■ 0 . 

By  (63)  this  is  only  possible  if  cot  w/m  - 0 , and  m > i , ie 

m - 2 or  2/3  . (69) 

The  form  of  the  condition  makes  it  very  likely  that  m - 2/(2M  + 1), 

M ■ 0,  1,  2,  ...  is  a complete  set  of  possible  values  of  m . However,  in 
subsection  4.1  the  possibility  that  m < } was  excluded  in  order  to  keep  the 
mapping  conceptually  simple,  so  the  present  argument  cannot  supply  a complete  set. 
In  fact  the  existence  of  the  tiro  possible  values  for  m , corresponding  by  (9) 
to  n ■ 3/2  and  5/2  in  (5),  and  the  non-existence  of  intermediate  values  is 
enough  for  practical  purposes.  It  is  suggested  in  section  7.3  that  larger  values 
of  n do  not  occur. 

4.4  Local  behaviour  of  velocity  and  pressure  on  the  wall 

We  now  consider  the  flow  on  the  wall  near  the  separation  line,  limiting 
the  discussion  to  the  two  values  of  m , given  by  (69),  for  which  a consistent 
solution  seems  to  be  possible.  The  velocity,  v , normal  to  the  separation  line, 
is  given  by  a form  of  (3): 


on  BA  and  DA,  where  v*  is  the  corresponding  velocity  component  in  the  Z^-plane 
and  is  given  by  (48)  and  (49),  with  (51)  for  BA  or  (53)  for  DA. 

1 /m 

On  BA,  C is  real  and  negative.  For  m • 2/3,  c'*“  - -ijc~/_|  , and  (20) 


gives 
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where  A^  + 0 and  Aj  tr*  und«t*rmin«d  constants.  Using  tha  ralAtions  abova 
to  express  C in  tarns  of  Z2  , And  writing  Zj  ■ y , w»  finally  find 


» - aj(i  *a?t  - t(- i)3/2  * 0<y2)) 

v ♦ 3(“  J)  + °(y>) 


for  n 
for  m 


2/3  1 


Thus  wa  saa  that,  for  n ■ 2 , as  tha  flow  approach* s tha  sapor at ion  lino, 
y - 0 , from  upstream,  tha  valocity  couponant  normal  to  tha  aaparation  lina  falls, 
and  its  gradiant  normal  to  tha  aaparation  lina,  v^  , bacomas  infinita  there. 

For  m • 2/3  no  such  bahaviour  arisas. 


On  tha  downatraaa  sida,  DA,  whara  u>  is  positiva,  (53)  with  (42)  givas: 


v* 

S 


maUg 

T55TT 


2sh-1 


TT 


wxu 


O + o(g*)) 


This  dominates  tha  ragular  part  of  v*  , given  by  (49).  Tha  modulua  of  tha 
mapping  is  given  by  (44)  and  so,  by  (70) 

* - nrrn;  <«  ♦ °<{2))  > 

or,  using  (28)  to  relate  g to  Zj  * y , we  have 

y ' -a.0?  ))*  <*  ♦ 0(y2/“»  . (74) 

As  would  be  expected  from  tha  proximity  of  the  sheet  to  the  wall  on  the  downstream 
sida,  this  expression  resembles  (64)  for  the  tangential  velocity  on  the  sheet. 

Tha  valocity  normal  to  tha  separation  line  on  its  downstream  side  is  zero  at  the 
lina  itself  and  directed  away  from  the  line  in  its  immediate  neighbourhood.  How- 
ever, it  is  so  small  that  tha  streamlines  on  the  surface  cross  the  conical  rays 
towards  tha  separation  lina,  since  v/U  < y/x  . Hence  the  velocity  normal  to  tha 
separation  line  decreases  along  a surface  streamline,  despite  the  appearance  of 
(74).  For  m ■ 2/3  or  2 the  velocity  in  the  cross-flow  plane  is  regular,  to 
tha  order  of  tha  present  analysis. 

We  now  examine  how  far  the  behaviour  of  the  pressure  gradient  near  the 
separation  line  reflects  that  of  the  gradient  of  the  velocity  component  normal  to 
the  separation  lina.  For  conical  flow,  the  pressure  gradient  along  the  separation 


line  is  t*ro.  Tbs  grtdisnt  normal  to  ths  ssparation  line,  from  the  upstream  to 
the  downstream  side,  is  3p/3y  . The  pressure  coefficient  is  given,  in  the  usual 
slander-body  approximation,  by 


2u  v2  ♦ w2 


+ const.  . 


On  the  wall,  s * 0 , where  w - 0 , 


- JPU2  ^ - - PCUUy  ♦ Wy) 


To  relate  uy  to  vy  , we  recall  that,  for  a conical  field, 


• - xF 


e-s 


and  so 


• • t&x  * y#y  ♦ s*#  - xu  + yv  + 


Differentiating  this  with  respect  to  y , 


xhy  ♦ v ► yvy  ♦ awy  , 


so  we  find,  for  a ■ 


- r\ 


and  the  pressure  gradient  becomes 


fc  - - ’%(*  - i)  ■ 


On  the  upstream  side,  where  v is  given  by  (73),  it  is  clear  that  the  pressure 
gradf  ent  behaves  in  the  same  way  as  the  velocity  gradient,  ie  there  ia  an  adverse 
gradient  which  becomes  infinitely  large  if  m ■ 2 ; while  if  m - 2/3  the 
gradient  is  not  locally  determined.  On  the  downstream  side,  with  v given  by 
(74),  the  pressure  gradient  becomes 


jR Ml 


(2m  + l)x 
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Since  the  pressure  itself  is  a function  of  y/x  in  a conical  flow,  it  is 
given  by 


P “ PD 


„2 

pU  m 

— ^ y 

(2m  + l)z 


(76) 


where  pD  is  the  pressure  on  the  separation  line.  Hence,  along  a surface  stream- 
line, on  which  as  we  have  seen  y/x  decreases  as  x and  y increase,  the 

pressure  falls  and  the  gradient  is  favourable.  With  a small  favourable  pressure 

gradient  and  a slowly  diverging  external  flow,  there  is  no  tendency  for  the 

boundary  layer  on  the  downstream  side  of  the  separation  line  to  separate.  This 
supports  the  concept  of  the  flow  advanced  in  section  2 and  sketched  in  Fig  1, 
regardless  of  the  value  of  m . 


A marked  difference  in  the  behaviour  of  the  limiting  streamlines  or  skin- 
friction  lines  between  the  upstream  and  downstream  sides  of  the  separation  line 
is  frequently  observed  in  surface  oil-flow  visualisations,  particularly  in 
approximately  conical  external  flows  with  highly-swept  separation  lines.  The 
limiting  streamlines  on  the  upstream  side  turn  towards  the  direction  of  the 
separation  line  very  late,  with  high  curvature  and  strong  convergence,  so  that 
their  behaviour  very  close  to  the  separation  line  is  often  masked  by  an 
accumulation  of  oil.  On  the  downstream  side,  the  limiting  streamlines  are  almost 
straight  and  parallel,  with  very  little  trend  towards  the  separation  line. 
Examples  are:-  primary  separation  on  a slender  circular  cone  at  incidence,  see 
for  instance  Nebbeling  and  Bannink' ; secondary  separation  beneath  the  primary 

g 

vortex  on  a slender  delta  wing,  see  for  instance  Lawford  ; and  both  primary  and 

9 

secondary  separations  on  an  upswept  rear  fuselage,  see  Peake  . Pictures  of 
glancing  interactions  (plane  shock  normal  to  the  surface  and  oblique  to  the 
stream)  by  Oskam,  Vas  and  Bogdonoff10  show  similar  differences  between  the 
upstream  and  downstream  sides  of  a region  of  streamline  coalescence,  although  the 
authors  do  not  regard  the  flow  as  separated. 

Although  the  behaviour  of  the  flow  on  the  downstream  side  of  the  separation 
line  is,  in  the  present  treatment,  little  affected  by  the  value  of  m , on  the 
upstream  side  the  value  of  m has  a marked  effect.  For  m - 2,  n ■ 3/2  , the 
existence  of  an  infinite  adverse  pressure  gradient  in  the  external  flow  makes  it 
probable,  though  not  certain,  that  the  boundary  layer,  whether  laminar  or  turbu- 
lent will  separate  before  it  reaches  the  postulated  separation  line.  This 
difficulty  does  not  arise  if  m • 2/3,  n - 5/2  . The  implications  of  this  for 
modelling  the  complete  flow  are  discussed  briefly  in  subsection  7.5. 
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5 EXTENSION  TO  NON-CON ICAL  SHEETS 

We  are  still  concerned  with  the  behaviour  of  a sheet  which  leaves  a plane 
wall  tangentially,  in  accordance  with  the  conclusion  of  section  2,  from  a separa- 
tion line  which  is  highly  swept,  ie  inclined  at  a small  angle  X to  the  undis- 
turbed  stream.  As  in  section  4 we  take  the  wall  to  be  the  plane  * ■ 0 , and 
choose  the  axis  of  x to  be  tangential  to  the  separation  line  at?  the  point 
considered,  with  x increasing  downstrsam.  We  again  study  the  configuration  in 
the  cross-flow  plane,  x - const  , and  use  slander-body  theory.  Near  the  separa- 
tion line  the  cross-section  of  the  sheet  is  again  of  the  form 

z - yyn  , n > I , (5) 


but  now  y and  n are  allowed  to  depend  on  x in  unspecified  ways.  It  will  be 
assumed  that  the  sheet  is  sufficiently  smooth  for  y and  n to  be  differentiable 
twice  with  respect  to  x . The  sequence  of  mappings  used  in  subsection  4.1  still 
applies  and,  by  (9),  the  variation  of  y and  n is  equivalent  to  some  equally 
smooth  variation  in  a and  m . It  is  again  assusMd  that  m > $ . 

The  boundary  conditions  for  this  system  of  axes  are  obtained  in  the 
Appendix,  where  the  normal  component  of  velocity  on  the  section  of  the  sheet  in 
the  cross-flow  plane  is  given  by  (A-5)  as 


v 

n 


sin  $ . 

6 const 


By  (8),  r 


a sinm(8/m)  and 


I 

r 3x 


so,  with  primes  denoting  derivatives, 

- m’  + m'  In  — + 0(82)  ; 

a m 


(77) 


so  that 

- m*  + m’  In  ~ + 0(62)^  r sin  ♦ . 

By  comparison  with  the  first  equation  of  subsection  4.2,  we  see  that  the  normal 
velocity  is  a multiple 


(78) 


of  the  conical  distribution,  plus  a contribution 
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vn  ■ - m'U  In  r sin  ♦ , (79) 

plus  a contribution  oJ  higher  order.  Recalling  that  the  procedure  followed  in 
•ubaection  4.2  to  obtain  the  tangential  components  from  the  normal  components  is 
linear,  we  see  that  the  tangential  velocities  in  the  non-conical  problem  can  be 
obtained  as  the  sum  of  (78)  times  those  found  for  conical  flow,  those  correspond- 
ing to  (79)  and  a contribution  of  higher  order.  The  last  will  be  neglected;  end 
an  analysis  like  that  of  subsection  4.2  will  be  pursued  for  the  normal  velocity 
given  by  (79).  The  tangential  velocities  so  found  will  be  used,  as  in  sub- 
section  4.3,  to  check  the  consistency  of  the  condition  of  pressure  continuity. 

The  same  geometrical  relations  as  were  used  in  subsection  4.2  enable  (79) 
to  be  rewritten  as 


vn  - - »'.U  .in”*'(^)  in  (|) 

■ - (f)<' * o",2/"»  • «*» 

in  place  of  (29).  The  geometrical  relations  (30)  to  (45)  still  apply,  so  we  can 
obtain  from  (80)  a set  of  equations  like  (47)  for  the  component,  w*  , of  the 
velocity  normal  to  the  real  axis  in  the  Z^-plane,  near  B and  D»  vis: 


on 

BA 

w*  - 

vn 

- 0 

* 

on 

BC 

w*  - 

K 

m'aU 

■ — nm  e>  - 

m 

(m+l)/m 

ln 

K1 

c(l  + 0(c1/n,  c ln  O) 

► 

on 

DC 

w*  - 

- V* 

ram'  aU  ^ 

2m+2  . , , . 

ln  6(1 

+ 0(62)) 

r. 

on 

DA 

w*  - 

V* 

- o . 

4 

n 


As  before,  we  regard  the  analytic  function  v*  - i w*  as  composed  of  a regular 

e 

part  and  a singular  part,  and  obtain  the  singular  behaviour  from  an  analytic 
function  whose  imaginary  part  behaves  like  v*  , as  given  by  (81). 

Starting  with  the  neighbourhood  of  B,  we  deal  first  with  the  possibility 
that  (m  + l)/m  is  integral,  ie.  in  the  admissible  range,  m ■ J or  1.  Consider 
the  function: 

f(;>  - - m c>2  . 


For  c real  and  negative,  </{f}  • 0 , for  5 real  and  positive, 
'ft | - in  c . Hnnea 


<»*  - i»*),  ~ - -=1&T  C(**l)/Ba*  - Id  C>: 
2 wink  | 


Adding  a regular  contribution,  aa  givan  by  (49),  loads  to 


v*  - iw* 


„n  m'aU  _(Brt‘l)/a/,  « wS2  . 

V - : rrsgi ' »*  - la  «>  ♦ •••  • 


noar  B,  for  m ■ } or  1,  where  Aq  is  roal  and  tha  tame  omitted  ara  of  highar 
ordar  than  tha  last  tarn  on  tha  right. 

Still  naar  B,  but  with  m and  & 1 , 1st  y • (m  + l)/m  and  consider 

«'>  - fe  - «'lu'  *» ')  • 

whara  is  defined  aa  when  K ■ pe*1^  . For  C real  and  negative, 

^ - w and  Jf\t\  ■ 0 . For  c real  and  positive,  /*f|  - In  C . Hence 


(v*  - iw*)s Sjgr  cosec  £ cC^l)M^-i*Ai  ln  c _ 


rr  cosec  — 1 . 


Adding  a regular  contribution,  as  before,  leads  to 


v*  - iw*  ■ 


00 

L*  _n  m’aU  ir  _(m+l)/m/-iir/m  , 

V - rsn- C0MC s'  (e  ln ' 

n-0  “*1 


- ir  cosec 


m)+  •’*  * 


. (83) 


near  B,  for  m > J and  =£  1 , where  A is  real  and  the  terms  omitted  are  of 

n 

higher  order  than  the  last  tern  on  the  right. 

Now,  near  D,  consider  the  function 


*k(“)  w ("  In  But)  ^ln  (“  In  But  / " * 


For  w real  and  negative, 
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(v*  - iw*)s  ~ (-  utsj)  (ln  (-  ttkz)  " sVr}  • (84) 

This  dominates  any  regular  contribution,  so  v*  - iv*  ~ (v*  - iw*)g  . 

The  next  step  is  to  obtain  the  tangential  velocities  on  the  sheet.  Consider 
first  the  upstream  side  BC,  near  B,  where  C is  small  and  positive.  (82)  and 
(83)  show  that  the  corresponding  component  in  the  Z^-plane  is 


v*  « v*  ■ An£n  ~ 
n»0 


(85) 


m’al)  C<m*l)/B  *2 

7m  2m+l  ln 
k1 


m'aU 


.(m+l)/m 


. tt  , 2m+l 
m sin  -■  k. 

IS  1 


From  this  vfc  follows  by  (54).  We  can  avoid  the  detailed  argument,  for  which 
the  various  values  of  m have  to  be  considered  separately,  by  observing  that. 
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although  the  singular  tarns  in  (85)  are  larger  than  those  in  the  previous  expres- 
sions (50)  to  (52) , they  are  still  of  higher  order  than  the  singular  terns  which 
appear  in  (dZ^/dZgj,  as  given  by  (32),  (34),  (36)  and  (38).  Consequently  the 
singular  terns  in  (85)  do  not  appear  in  the  expressions  for  vt  , which  are 
identical  with  (55)  to  (58).  Again,  the  dominant  terms  in  v near  B are  the 
same  as  they  would  be  if  the  sheet  were  a streamline  of  the  cross-flow. 

On  the  downstream  side  DC,  near  D,  on  the  other  hand,  the  tangential  com- 
ponent is  not  the  same  as  in  the  case  of  conical  flow.  The  corresponding  velocity 
component  in  the  Z ..-plane,  v*  , is  -v*  - -v*  , as  given  by  the  real  part  of 
(84)  for  w real  and  negative: 


mm'aU 

vt  " it  (2m  + 1)« 


s2t  '‘(ln  6 ’ ^tt)(1  * 0<{2>)  • 


Hence,  by  (43)  and  (54) 


sttt  *“(*» 8 - srnr)(1  + °<s2»  • 


or,  by  (45), 


m*Uo  f,  a m \ ,,  . 2/mXv 

ssmy  i - a-rr)  (1  ♦ 0(0  » • 


(86) 


This  is  of  lower  order  than  the  corresponding  result  (59)  for  the  conical  sheet. 

We  can  now  write  the  tangential  velocity  components  on  the  sheet  near  B 
and  near  D for  the  complete  non-conical  problem  by  combining  (78)  times  the 
expressions  for  conical  flow  with  those  just  obtained  for  the  special  problem 
defined  by  (79).  Near  B,  both  expressions  are  the  same,  so  the  result  has  still 
the  form  found  for  the  conical  flow: 


rfl  + Aj'o  + (m  + D(f)1/mf (o)  + o(o,/ffif(a))j  , 


vt'BC  * A0 


where  Aq  # 0 and  Aj’  are  undetermined  constants  (not  derivatives!)  and 


f(o) 


In  a 


,cot  (=) 


for  m - 4 or  1 


for  m > | and  ^ 1 


(62) 


(63) 


> ■ "if  i u -u-i 
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Near  D,  by  (86)  and  (64) , we  have 


, Ua  /m'  a , a'  2(m  + l)m'\  ,,  . n,_2,  aXK  ,M* 

'tine  2nT+T  ("in"  1 a 1 STT  ) (1  0(0  » • <*” 


The  condition  of  continuity  of  pressure  across  the  sheet  is  given  in  the 
Appendix  as 


Mi 

3x 


!8  const 


- iv/| s. 
tl  3x 


cos  4 


'6  const 


■>)■ 


(A-ll) 


From  (62)  and  (87)  we  can  write 


Avt  - a/.  ♦ C,0.  Cm  + 1)  ln  f * •••  • ««> 


where  Cj  is  another  undetermined  constant  and  the  omitted  terms  are  of  higher 
order  than  the  smallest  term  included.  Integrating  (88)  along  the  sheet  gives 


a*  - r 


(f  o ,_v<m+0/tt  \ n_»„2 

« * V ♦ “(!)  t(a))  ' 2S?sV 


7y  in  - + ... 


(89) 


where  r is  again  the  total  circulation  of  the  sheet,  is  a constant  related 

to  Cj  and  the  omitted  terms  are  of  higher  order.  This  equation  expresses  A* 
in  terms  of  a and  a number  of  quantities,  A^,  C^,  m and  a , which  are  func- 
tions of  x . It  is  therefore  helpful  to  write: 


Mil 

'0  const 


Mil  + Mi 

3a 


3x 


<o  const 


3a 
dx 

x const  '6  const 


Introducing  this  into  (A- 11),  we  have 


Mi 

3x 


a const 


6 const 


cos  ♦ - T- 


3a 

3x 


m 


9 const 


U 


(90) 


Now  a " a(0/m)m(l  + 0(0^))  , so 


3a  j 
3x  L 

'8  const 


m a \ a / 
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while,  by  (61)  and  (77), 


M|  cos  ♦ 

1 0 const 

Bence 


■ (r"  lnf-t(T‘m’)<') 


(1  + 0 (o2/o)) 


3r  , 3o 

JZ  COS  ♦ - T- 

8 const  8 const 


. 0(o(m*2)/"‘  In  c>  . 


(91) 


Bence,  by  (89)  and  (91),  we  can  write  (90)  as 


a / 2 (m+l)/m 

2 tt-  IT  + AAo  + A’C2o  + A^ma^-J 


3x 


]T  ♦ Ajc 


Ura'o  a . 

f(o;  " 2m(2m  + 1)  iU  a *** 


- ^vc ! DC' 2 ” ^vJkP2  + 0(o(n+2)/°  ln  °)  » (92) 


where  the  derivative  is  for  constant  a and  the  velocities  are  given  by  (87)  and 
(62). 

k2 


Now,  by  (62),  (vtl rc'  contains  a term 


1/m 


|A  « / IU 

2A0Z(m  + l)(SLj  f (o)  . 


(93) 


If  m > 1 and  ¥*  1 , the  exponent  of  o in  (93)  is  fractional  and 

f(o)  - cot  (ir/m)  . No  other  term  in  (92)  involves  a to  a fractional  power  as 

small  as  1/m  , so  (93)  roust  be  zero.  Since  A^  =£  0,  f (a)  - cot  (ir/m)  » 0 , 

and  so  m • 2/3  or  2 . Hence  m cannot  take  values  which  are  close  to  } or  1 ; 

but  m is  continuous,  so  m can  only  take  the  values  | or  1 if  m*  ■ 0 . 

2 2 

If  m - | , the  term  (93)  is  3A^  (o/a)  In  a . The  only  other  terms  in  (92)  of 

order  o^  In  o have  m*  as  a factor,  so  there  is  nothing  to  balance  (93).  The 

term  cannot  vanish  and  so  m ■ i is  excluded.  If  m - 1 , (93)  is 
2 

4A^  (o/a)  In  a and  there  is  no  other  term  of  this  order  in  (92).  Hence  m - 1 
is  also  excluded. 

The  only  possible  values  for  m are  therefore  2/3  and  2,  just  as  in  the 
case  of  the  conical  vortex  sheet. 

Just  as  the  dominant  terms  in  the  tangential  velocity  on  the  upstream  side 
of  the  sheet  are  not  affected  by  allowing  the  sheet  to  be  non-conical,  wo  can 
see  that  the  dominant  terms  in  the velocity  and  pressure  on  the  wall  upstream  of 
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the  separation  line  are  also  unchanged.  There  are  again  adverse  gradients  of 
cross-flow  velocity  and  of  pressure , becoming  infinite  at  the  separation  line,  if 
m ■ 2 j but  no  locally  determined  gradient  if  m - 2/3  . On  the  wall  on  the 
downstream  side,  a behaviour  similar  to  that  of  the  conical  flow  will  be  found 
again,  since  the  additional,  lower-order  term  in  the  velocity  field,  (84), 
vanishes  for  m'  - 0 . The  details  will  be  more  complicated,  owing  to  the  extra 
freedom  available  in  non-conical  flow,  and  will  not  be  explored  here. 

6 EXTENSION  TO  CURVED  WALLS  IN  NON-CONICAL  FLOW 

In  this  section  the  restriction  that  the  wall  should  be  plane  will  be 
replaced  by  the  weaker  restriction  that  it  should  be  smoothly  curved.  The  argu- 
ment is  arranged  so  as  to  exploit  the  results  obtained  already  and  avoid  repeti- 
tion of  the  analysis. 

The  framework  of  slender-body  theory  is  retained,  so  the  curved  wall  is  the 
surface  of  a slender  body  at  a small  angle  of  incidence  to  the  uniform  stream. 

The  separation  line  on  it  is  assumed  to  be  smooth  and  highly  swept,  so  that  there 
is  a tangent  to  it  at  the  point  considered  and  this  tangent  is  inclined  at  a small 
angle  to  the  uniform  stream.  We  choose  this  tangent  as  the  axis  of  x , with  x 
increasing  downstream.  The  body  is  also  supposed  to  be  smooth  at  this  point,  so 
it  has  a tangent  plane,  and  we  choose  this  as  the  plane  z.  - 0 . with  z . 
increasing  away  from  the  body.  The  y j-axis  is  chosen  so  that  (x,yJtZj)  form  a 
right-handed  system.  Then,  according  to  slender-body  theory,  the  flow  in  the 
plane  x - constant  may  be  described  in  terms  of  analytic  functions  of 
Z,  - y,  ♦ i.,  . 

The  argument  is  again  a local  one,  so  that  the  sheet  and  body  surface  can  be 
replaced  by  other  curves  which  have  the  same  local  behaviour  near  the  point  con- 
sidered. The  cross-section  of  the  smoothly-curved  wall  can  be  replaced  locally 
by  a parabola: 

z j ■ - hy*  , h - h(x)  . (94) 

Then  the  conformal  mapping 

Zj  - Z2  - ihZ*  (95) 

maps  the  region  above  the  parabola  (94)  in  the  Zj -plane  on  to  the  upper  half  of 
the  Z ^ -plane . We  regard  this  as  the  region  occupied  by  the  fluid,  so  that  the 
body  is  locally  convex  if  h > 0 . The  Zj  and  Zj-p lanes  are  sketched  in  Fig  5. 
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On  the  same  basis,  the  sheet  is  replaced  by  the  same  curve  in  the  Z. -plane  that 
was  used  in  sections  4 ana  5 to  represent  the  sheet  in  that  planet  ie 

Z2  " r2ei02  » r2  " a * 0 < 02  « «ec  , (96) 

as  in  equation  (8).  The  same  sequence  of  transformations  as  bafore  maps  the 
region  exterior  to  the  body  and  sheet  on  to  the  upper  half  of  the  Z^-plane,  but 
we  shall  not  need  to  make  explicit  use  of  these. 

We  must  first  check  that  a general  choice  of  the  parameters  a and  m 
(and  of  their  variation  with  x)  provides  the  desired  degree  of  generality  of 
sheet  shape  in  the  cross-flow  plane.  If  the  distance  of  the  sheet  from  the  body, 
measured  parallel  to  the  Zj-axis,  is  denoted  by  f(y{)  , then  by  equations  (94)  to 
(96): 


yj  ♦ if(yj)  - Zj  ♦ ihyZ  - Z2  - ihZZ  ♦ ibyZ 


ifl2  ..  2 2i02  , • » 2 

- r2e  *■  - lhr^e  1 + ihyj  . 


Hence,  for  small  values  of  82  , 


yj  » r2  cos  82  + hr2  sin  202  ~ aj 


and 


f (y,)  • r,  sin  8,  - hr,  cos  28_  + 


© 

"’i  ~ -(tJ 


m+1 


(m+l)/m 


This  corresponds  to  the  general  form  expressed  by  equation  (5)  if 

.m 


m 


1 

n - 1 


and 


a 


-tf) 


just  as  in  (9).  The  form  given  by  (96)  is  therefore  sufficiently  general. 

The  first  stage  of  the  argument  is  to  relate  the  normal  velocity  on  the 
sheet  in  the  Zj-plane  which  is  required  to  make  it  the  section  of  a stream 
surface  of  the  three-dimensional  flow  to  the  normal  velocity  imposed  on  the 
corresponding  curve  in  the  Z2~plane  in  the  previous  sections.  At  the  same  time 
we  treat  the  normal  velocity  on  the  wall.  The  stream  surface  boundary  condition 


Is  expressed  by  equation  (A-4 ) , in  which  the  parameter  held  constant  in  forming 
the  x-derivatives  may  be  chosen  at  will. 


On  the  wall  (which  equals  y,  , by  (95))  is  held  constant,  so  that 

(A-4)  becomss 


from  which 
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where  vft,  and  are  the  normal  velocities  in  the  Z,  and  Z ^ planes. 

a ^ ■ y2  “ y,  on  the  wall  and  z,  is  given  by  (94),  so 


n2 

D 


“ h*yo 


On  the  sheet  82  is  held  constant,  so  that  (A-4)  gives 


v»2  _ vnl  dt,l 


9*. 


U U do2  3x 


ay,  ay, 


do2  8x 


dz. 


do. 


Inserting  (96)  into  (95)  gives 


y,  - r2  cos  8j  ♦ hr 2 sin  202 


!j  - r2  sin  02  - hr^  cos  2©2  , 


from  which 


3r„ 


3x 


(cos  ©2  + 2hr2  sin  202>  + h'r*  sin  202 


3z 


3r„ 


r-r  - (sin  8,  - 2hr,  cos  20.)  -r-r-  - h’ri?  cos  20 


(97 


(98. 


(99] 
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Also,  remembering  that  cos  *2  - dr2/d o2  and  sin  - e2d02/do2  , and  that 
♦2  * 82  ♦ 42  , whara  ia  tha  Inclination  of  tha  aaction  of  tha  shaat  to  tha 
raal  aids,  va  hava 


35-  - coa  *2  ♦ 2hr2  ain  («2  ♦ *2> 


35!  - ain  *2  - 2hr2  coa  («2  > *2>  . 


(100) 


Introducing  (99)  and  (100)  into  tha  boundary  condition  above  yield* 


- - ^1  ♦ 4hr2  ain  e2  ♦ 4h2r2)  ain  ^ 
- h'r2(coa  (fl2  - $2)  + 2hr2  ain  . 


(101) 


To  find  tha  dominant  terms  in  (101),  all  tha  quantitias  need  to  bo  expressed  in 
tarms  of  one,  say  02  . By  (96) 

^ *2\* 


d02  2 /dr2\ 

d6^  r2  l d0^J 


. m-I  2 

- a am  •—  , 

ID 


so  that  ai  >2  <*  r2d02/do2  - ain  (02/m)  , 


♦2  - r 


a 

-*) 


002) 


(103) 


dr9/3x  ven  > /)  with  the  omitted  suffix  2 restored  and  r0  itself  is 

^ 2 2 2 /m  ^ 

given  by  (96).  Neglecting  92  ■-  tn  (a^/a)  in  comparison  with  unity,  as  in  the 

previous  analysis  (of.  (29)  and  (80)),  ve  can  write  (101)  as 


■ &) 


In  ~ + a’  - am'l  1 + 4h« 
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1 * 2h*hfl 
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From  this  we  can  pick  out  tha  dominant  tacma  in  (101 ) t 


•in  42 


004) 


Hhich  of  thaaa  two  dominataa  dapanda  on  tha  valua  of  a,  to  both  ara  ratainad, 
using  a small  extans ion  of  tha  tilda  notation. 

Tha  normal  velocitias  spacifiad  by  (97)  and  (104)  ara  tha  sum  of  two 
contributions : 

v_2  3*2 

(a)  vn2  ■ 0 on  tha  wall  and  sin  42  on  ^ ■hast,  and 

V V 

(b)  • -h'y^  on  the  wall  and  - -h'o^  on  tha  sheet. 

Tha  linearity  of  the  relation  between  the  normal  and  tangential  velocity  distribu- 
tions can  be  exploited,  as  it  was  in  section  5,  so  that  the  complete  tangential 
velocities  can  be  obtained  as  the  sum  of  those  corresponding  to  distributions  (a) 
and  (b) , Now  distribution  (a)  is  precisely  that  treated  in  sections  4 and  5, 
arising  from  the  sheet  on  the  plane  wall.  For  distribution  (b) , we  guess  a 
complex  potential 

^ - ih*UZ^3  (105) 

and  verify  that  it  gives  the  right  behaviour  on  the  wall  and  sheet.  For  the 
wall,  it  is  enough  to  differentiate  and  set  Z2  - y2  , real: 


This  verifies  the  behaviour  on  the  wall.  For  the  contribution  from  on  the 
sheet,  it  is  convenient  to  write 


dW^  dW^dZj 

vt2  " lvn2  “ do2  " dZ2  do2  * 

Now  dZ2/da2  ■ e^2  - e*-(m+^®2/m  f by  (102).  Bence 

vt2  - ivn2  - ih’Ur2ei(3B+1)02/,n  , (106) 


so  thst 
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and  tha  behaviour  on  the  sheet  is  verified,  to  the  usual  approximation.  Bence, 
to  leading  order,  the  tangential  velocities  in  the  ate  those  found  in 
the  previous  section,  augmented  by  those  derived  from  (105). 


The  final  step  is  to  consider  the  compatibility  of  these  tangential 
velocities  with  the  condition  of  continuity  of  pressure  across  the  sheet.  This 
condition  can  be  written  in  the  form  (A- 10) , and  it  is  again  convenient  to  choose 
the  parameter  n to  be  82  • Then: 
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1 m Aw  1 
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(107) 


where  the  subscripts  '1'  indicate  quantities  in  the  Zj -plane.  In  section  5 an 
expression,  (89),  was  found  for  4#  in  terms  of  a . Sin-.*  the  extr*.  velocity 
(106)  is  continuous  across  the  sheet,  this  expression  still  holds,  to  the  order 
considered,  with  o replaced  by  O2  • Ifc  i*  therefore  convenient  to  consider, 
as  before,  its  derivative  for  constant  c?2  , and  write 
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The  velocity  componerts  in  the  Zj  and  Z^ “planes  are  related  by 


(108) 


Using  equations  (99)  and  (100)  and  the  relation 


1 + 4hr2  sin  82  + 


4h2r2 


(109) 
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which  holds  on  the  those,  vs  find 


010) 


Introducing  (110)  into  (108)  and  r seal ling  (91),  vs  havo 

.2 


ir|o  • ^2(357)  H<*in  <"2  - ♦ 
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The  additional  tangential  component  of  velocity  deriving  from  is  given  by 
(106),  so  that 


Vtm2  Vtm  ^,  2 £n  + 1 „ 
_ - — - h r2  am  e2 


and 


&vt2  " Avt  * 


where  v and  Av  are  just  as  they  were  in  section  5.  Using  (103)  and  (109) 

til  t 

again,  we  find  the  condition  of  pressure  continuity  becomes 

- Avt(l  + 4hr2  sin  02  + 4h2r2)^-  ♦ h'r2  sin  «2 


+ h'r?  sin  8*  + 2hh*rJ  cos 


l2  m 2 2 

The  corresponding  equation  in  section  5 is  (92),  which  can  be  written 
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013) 


• ! /m 
It  was  shown  that  this  contains  an  unbalanced  term  of  order  a if  m > | 

and  m # 1 or  o^11  In  o if  m - J or  1.  Recalling  that  Avt  and  v^  are 

both  of  order  unity,  we  can  see  that  the  orders  of  the  extra  terms  which  appear 

in  (112)  and  not  in  (113)  are  too  high  for  them  to  balance  this  crucial  term. 

It  must  therefore  vanish,  and  the  argument  in  section  5 which  follows 

equation  (92)  again  holds,  showing  that  the  only  possible  values  of  m in  the 

range  $ 4 m < « are  still  2/3  and  2. 
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The  extra  freedom  afforded  by  the  smoothly  curved  wall  is  not  therefore 
sufficient  to  allow  a different  behaviour  of  the  sheet  close  to  the  separation 
line.  Further,  since  the  additional  velocity  field  defined  by  (105)  and  (95)  is 
regular,  the  singularity  in  the  pressure  gradient  on  the  wall  found  previously 
on  the  upstream  side  of  the  separation  line  for  m - 2 (n  ■ 3/2)  will  still 
arise. 

The  shape  of  the  sheet  is  described  most  simply  in  terms  of  its  distance 
f(yj)  from  the  wall,  y j being  the  same,  to  the  first  approximation,  as  the 
distance  along  the  wall.  Wo  have 

fCy,)  ~ . 


where  n ■ 1 + I/m  is  3/2  for  m » 2 and  5/2  for  m - 2/3  . In  terms  of  the 
Cartesian  coordinates  y^  and  Zj  measured  parallel  and  perpendicular  to  the 
tangent  to  the  body  cross-section  (see  Fig  5),  the  shape  is  given  by  equations 
(95)  and  (96)  as 


z,  + hy 


for  m ■ 2 


for  m ■ 2/3  . 


7 RESULTS  AND  CONJECTURES 
7. 1 Main  results 

It  has  been  shown  that,  under  the  assumptions  of  slender-body  theory,  the 
shape  of  a vortex  sheet  which  contains  circulation  being  shed  from  a highly-swept 
separation  line  on  a smoothly  curved  wall  in  steady  flow  is  restricted  in  form. 
If,  in  the  vicinity  of  the  separation  line,  the  distance  of  the  sheet  from  the 
wall  increases  like  y11  , where  y is  the  distance  from  the  separation  line, 
then  the  only  possible  values  of  n not  exceeding  3 are  3/2  and  5/2.  If 
n - 3/2  , there  is  an  adverse  pressure  gradient  on  the  wall  upstream  of  the 
separation  line  which  becomes  infinite  at  the  separation  line  itself.  If 
n - 5/2  , the  pressure  gradient  at  the  separation  line  remains  finite.  In  both 
cases  the  strength  of  the  vortex  sheet  (its  circulation  per  unit  length)  varies 
linearly  with  the  arc  length  along  the  sheet,  close  to  the  separation  line. 
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For  a conical  ahaac  on  a plana  wall,  tha  inviacid  flow  on  tha  downstream 
side  of  tha  aaparation  lina  it  favourabla  to  tha  growth  of  an  attached  boundary 
layer,  to  that  it  it  only  tha  boundary  layer  on  tha  upstream  tide  which  teparates. 

Theta  retultt  are  expected  to  be  helpful  in  numerical  treatment!  of  vortex 
sheet  models  of  separated  flow  and  in  understanding  the  mechanics  of  flow 
separation  in  three  dimensions. 

7.2  Equivalent  results  for  unsteady  plane  flow 

There  is  an  exact  analogy  between  slander-body  theory  for  steady  flow  and 
the  theory  of  the  time-dependent,  two-dimensional  flow  of  an  inviscid  incompress- 
ible fluid.  The  time,  t , in  the  plane  flow  is  related  to  the  downstream  dis- 
tance, x , in  the  steady  three-dimensional  flow  of  undisturbed  speed  U by 


x 

U 


(114) 


;v 


fe  I 


The  plane  flow  takes  place  in  the  cross-flow  plane.  Variations  in  the  cross- 
section  of  the  slender  body  correspond  to  changes  with  time  in  the  shape  of  the 
body  in  the  plane  flow  and  variations  in  its  inclination  to  the  undisturbed 
stream  correspond  to  changes  in  the  relative  velocity  between  the  body  and  the 
fluid  in  the  plane  flow.  The  term  3*/3t  which  arises  in  the  time-dependent 
form  of  Bernoulli's  equation  corresponds  to  the  contribution  from  the  streamwise 
perturbation  velocity  u - 3$/3x  to  the  slender  body  form  (75)  of  the  three- 
dimensional  ' Bernoulli  relation.  It  is  straightforward  to  verify  that  the 
boundary  conditions  on  the  vortex  sheets  correspond  through  equation  (114). 

The  results  of  the  present  work  therefore  hold  also  for  the  behaviour  of  a 
vortex  sheet  representing  an  'open'  type  of  separated  flow  in  plane,  incompress- 
ible time-dependent  flow,  whether  the  time  variations  are  in  the  imposed 
velocity,  the  body  shape  or  the  position  of  the  separation  point.  In  particular, 
the  same  two  exponents  3/2  and  5/2  are  the  only  ones  (not  greater  than  3)  which 
can  occur  in  flows  of  this  kind,  and  the  smaller  is  inevitably  associated  with 
an  infinite  adverse  pressure  gradient. 

7.3  Relation  to  previous  results  for  plane  flow 

A considerable  amount  of  work,  dating  back  to  Helmholtz  and  Kirchhoff, 
has  been  done  on  steady  two-dimensional  flows  past  a body  from  which  spring  a 
pair  of  constant-pressure  streamlines.  Since,  in  this  context,  constant  pressure 
implies  constant  speed,  hodograph  methods  can  be  used  and,  with  their  help,  com- 
plete flowfields  can  be  described  analytically.  A summary  account,  with 
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references  to  the  original  papers,  is  givan  by  Thwaites 1 1 . Tha  cons tan t-prasaura 
atraaalinaa  laava  tha  body  aurface  tangantially  and,  in  soma  Sanaa,  anclosa  a 
finita  or  infinita  region  downstream  of  tha  body.  Tha  fluid  inside  this  region 
must  be  at  rest  at  the  separation  point  itself,  since  it  lias  between  the 
constant-pressure  streamline  and  tha  downstream  side  of  tha  body,  which  meat  at 
a cusp.  Since  tha  static  pressure  is  continuous,  the  constant  pressure  on  tha 
separation  streamlines  must  be  tha  stagnation  or  total  pressure  in  tha  downstream 
region.  Tha  total  pressure  is  therefore  different  in  tha  downstream  region  from 
its  value  in  tha  main  flow  and  the  separated  flow  cannot  be  of  the  'open*  type 
discussed  in  the  present  Report.  There  are  two  conventional  interpretations  of 
the  downstream  region:  either  the  fluid  there  is  the  same  as  the  mainstream 
fluid  and  it  is  at  rest  (a  deadwatcr  region),  or,  more  plausibly,  it  is  occupied 
by  a fluid  whose  density  is  negligible  compared  to  that  of  the  mainstream  fluid 
(a  cavity  in  a liquid,  filled  by  its  vapour  or  another  gas).  Neither  of  these 
interpretations  seems  to  have  any  aerodynamic  significance  in  itself;  but,  since 
the  constant-pressure  streamlines  are  particular  cases  of  vortex  sheets  springing 
from  a smooth  body,  their  behaviour  near  the  separation  point  invites  comparison 
with  the  present  results. 


In  the  vicinity  of  the  separation  point,  equations  (168)  and  (169)  of 
Ref  11  describe  the  behaviour  of  the  curvature,  k , of  the  streamline  and  the 
upstream  pressure  gradient  on  the  wall.  In  condensed  form: 
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where  T is  a constant,  s is  the  distance  along  the  body  surface  in  the  direc- 
tion of  the  upstream  flow  and  4 is  the  velocity  potential  on  the  sheet  or  body, 
taking  the  value  *2  at  the  separation  point.  We  can  see  that  if  T > 0 the 
curvature  and  pressure  gradient  behave  in  the  same  way  as  the  present  results 
predict  for  n ■ 3/2  . If  T < 0 , the  solution  is  meaningless,  as  the  constant- 
pressure  streamline  lies  initially  within  the  body.  If  T « 0 , both  the  curva- 
ture and  pressure  gradient  vanish.  The  present  results  are  that  the  curvature 
and  pressure  gradient  are  finite  for  n - 5/2  ; the  curvature  is  that  of  the  body 
and  the  pressure  gradient  is  not  locally  determined.  Hence  the  previous  results 
for  the  neighbourhood  of  the  separation  point  reduce  almost,  but  not  quite,  to  a 
special  case  of  the  present  analysis.  The  difference  concerns  the  curvature  of 
the  sheet  in  the  case  for  which  it  is  finite:  the  constant-pressure  streamline 
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has  zero  curvature,  the  'open'  type  of  vortex  sheet  has  the  same  curvature  as  the 
body.  It  is  plausible  that  such  a difference  should  arise,  since  the  flow 
between  the  sheet  and  the  body,  on  the  downstream  side  of  the  separation  point, 
does  not  enter  the  previous  analysis. 

The  previous  results  also  suggest  the  role  played  by  the  solutions  for 
n - 3/2  and  5/2.  To  see  this,  we  confine  our  attention  to  the  so-called 
Kirchhoff  flow  past  a circular  cylinder,  the  flow  for  which  the  pressure  on  the 
constant-pressure  streamline  is  the  same  as  the  pressure  in  the  undisturbed 
stream,  as  illustrated  in  Fig  IV. 14  of  Ref  11.  There  is  a one-parameter  family 
of  these  flows,  corresponding  to  a variation  in  the  position  of  the  separation 
points  on  the  cylinder.  If  the  position  of  separation  is  specified  by  its 
angular  displacement  6 from  the  front  stagnation  point,  then  for  6 > 0g  « 56°  , 
it  is  found  that  T > 0 . As  6 falls  towards  0,,  , T falls,  and  vanishes  when 
0 * 0g  . Separation  which  occurs  with  T - 0 is  called  'smooth'.  Suppose,  now, 
that  we  find  a vortex  sheet  representation  of  the  separated  flow  past  a slender 
cone  (of  general  cross-section)  with  an  assumed  inviscid  separation  line  along 
a generator  chosen  at  random;  then  it  3eems  likely  that  the  behaviour  near  the 
separation  line  will  correspond  to  n - 3/2  . The  behaviour  corresponding  to 
n - 5/2  is  not  likely  to  arise  for  more  than  a finite  number  of  positions  of  the 
separation  line,  though  such  positions  may  well  be  significant  in  relation  to  the 
real  flow,  as  discussed  in  subsection  7.5.  Vortex  sheets  which  lie  even  closer 
to  the  surface,  which  might  arise  from  values  of  n greater  than  3,  seem 
unlikely  to  occur  at  all. 

7. A Range  of  validity  of  present  results 

The  detailed  results  of  the  present  work  have  been  obtained  within  the 
assumptions  of  slender-body  theory.  It  is  relevant  to  ask  whether  they  hold  more 
widely. 

On  the  one  hand,  the  proof  in  section  2 that  a vortex  sheet  representing 
separation  must  be  tangential  to  the  wall  depends  crucially  on  the  equality  of 
the  entropy  on  the  two  sides  of  the  sheet.  This  clearly  sets  a limit  to  the 
extent  to  which  the  later  results  about  the  flow  near  the  separation  line  can  be 
generalized.  On  the  other  hand,  since  the  arguments  used  in  obtaining  these 
results  are  local  ones,  it  seems  unlikely  that  the  global  limitations  required 
for  the  application  of  slender-body  theory  are  actually  necessary.  This  impres- 
sion is  confirmed  by  the  applicability  of  the  present  results  to  plane  flow,  a 
two-dimensional  body  being  completely  ' non-s lender ' . Further,  the  applicability 
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of  the  present  results  to  unsteady  plane  flaw  suggests  that  the  initial  limita- 
tion to  steady  flow  was  unnecessary. 

It  seems  that  it  might  well  be  possible  to  extend  the  present  results  to 
unsteady  incompressible  flow  past  general  bodies,  and,  perhaps,  to  compressible 
flows  in  small-disturbance  approximations  in  which  entropy  changes  are  negligible. 

7.5  Prediction  of  separation  line 

Studies  of  inviscid  flew,  like  the  present,  are  meant  to  lead  to  useful 
models  of  real  flows.  As  suggested  in  section  1,  vortex  sheet  models  of  flows 
separating  at  highly-swept  separation  lines  are  expected  to  be  adequate  if  the 
position  of  the  separation  line  can  be  determined.  In  Ref  3 it  was  suggested 
that  this  should  be  attacked  by  an  iteration  between  a boundary-layer  calculation 
which  uses  a given  external  velocity  field  to  predict  the  location  of  separation 
and  a calculation  of  an  inviscid  model  of  the  separated  flow  which  uses  a given 
separation  line  to  predict  the  external  field.  For  turbulent  boundary  layers 
this  still  appears  to  be  the  way  to  proceed,  but  for  laminar  flows  there  may  be 
an  alternative  procedure. 

The  resemblance  between  the  present  results  for  separation  from  highly- 

swept  separation  lines  and  earlier  results  for  the  constant-pressure  streamline 

12 

model  of  two-dimensional  separation  suggests  that  an  analogue  of  Sychev's  work 
for  plane  flows  may  exist  for  highly-swept  separation  lines.  Sychev  seeks  a 
local  solution,  valid  near  the  separation  point,,  to  the  Navier-Stokes  equations 
for  the  steady  flow  of  an  incoaq>ressible,  viscous  fluid  past  a smooth  body  in 
two  dimensions.  He  finds  the  structure  of  such  a solution  in  the  form  of  an 
asymptotic  expansion  in  fractional  powers  of  the  inverse  of  the  Reynolds  number. 

As  the  Reynolds  number  tends  to  infinity,  the  solution  tends  towards  an  inviscid 
flow  with  a constant-pressure  streamline  which  separates  smoothly  from  the  wall. 
For  finite  Reynolds  number,  the  point  of  separation  is  further  downstream,  by  a 
distance  proportional  to  Re  . The  constant  of  proportionality  remains  to 

be  determined  by  a numerical  calculation,  which  is  required,  in  principle,  to 
confirm  that  a solution  with  this  structure  actually  exists. 

If  such  an  analogue  of  Sychev's  work  does  exist,  and  preliminary  work  by 
N.  Riley  at  the  University  of  East  Anglia  indicates  that  it  may  do  so  for  conical 
external  flows,  then  the  outcome  of  the  iteration  proposed  above  is  predictable: 
since  laminar  boundary  layer  theory  represents  the  limit  of  infinite  Reynolds 
number,  it  is  smooth  separation  which  will  emerge.  This  would  be  the  best  that 
could  be  achieved  with  a bound ary- layer  calculation  and  it  could  be  achieved 
without  it,  by  examining  the  behaviour  of  the  inviscid  model.  Furthermore,  if 
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the  numerical  solution  to  the  analogue  of  Sychev's  problem  vere  known,  the  lead- 
ing term  (of  order  Re  ^ sS)  in  the  correction  to  this  infinite  Reynolds  nuaber 
solution  could  be  found,  and  the  inviacid  model  corresponding  to  this  corrected 
position  of  separation  could  be  calculated  for  large,  but  finite,  Reynolds 
nunbers. 

This  is  all  speculative  at  present,  but  already  has  a bearing  on  the  con- 
struction of  an  inviscid  model  of  the  flow  separating  from  highly-swept  separa- 
tion lines  on  smooth  surfaces.  If  such  a model  is  intended  eventually  to  form 
part  of  a method  which  predicts  the  position  of  the  separation  line, 

(a)  it  must  represent  the  shed  vorticity  near  the  separation  line  in  a 
continuous  form,  ie  a vortex  sheet  representation  is  the  simplest  likely 
to  be  useful, 

(b)  the  curvature  of  the  sheet  must  be  allowed  to  tend  to  infinity  like  the 
inverse  square  root  of  the  distance  from  the  separation  line,  and 

(c)  the  coefficient  of  this  singularity  should  be  determined  explicitly  by  the 
calculation  procedure. 

Finally,  since  the  constant-pressure  streamline  model  has  figured  so 
prominently  in  this  section,  its  role  should  perhaps  be  clarified.  The  resem- 
blances between  the  mathematical  properties  of  the  vortex  sheet  model  of  open 
separation  (whether  in  steady  flow  from  a highly-swept  separation  line  or  in 
unsteady  two-dimensional  flow)  and  the  constant-pressure  streamline  model  of 
steady  two-dimensional  separated  or  cavity  flow  have  been  noted  and  partially 
exploited.  It  is  not  suggested  that  the  constant-pressure  streamline  model  has 
the  same  relevance  to  closed  separation  in  two-dimensional  steady  flow  as  the 
vortex  sheet  model  is  believed  to  have  to  open  separation.  Indeed,  a more 
satisfactory  approach  to  two-dimmnsional  separation  would  be  to  recognize  its 
essentially  unsteady  character  and  treat  it  as  a time-dependent  problem  involving 
vortex  sheets. 
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In  this  Appendix  the  boundary  conditions  on  the  body  and  on  the  steady 
vortex  sheet  are  expressed  in  terms  of  quantities  in  a cross-flow  plane,  under 
the  usual  assumption  of  small  disturbances.  The  expressions  obtained  are  equiva- 
lent  to  those  given  by  Smith  and  Clarke  ; but  their  derivation  is  a little  more 
general,  in  that  the  axis  system  is  no  longer  required  to  be  exactly  aligned  with 
the  undisturbed  stream,  and  a little  simpler,  in  that  vector  algebra  is  avoided. 

The  boundary  conditions  are  that  the  body  and  the  vortex  sheet  are  stream 
surfaces  of  the  three-dimensional  flow  and  that  the  pressure  is  continuous  across 
the  vortex  sheet.  Since  the  sheet  is  assumed  to  be  embedded  in  an  irrotational 
flow,  the  continuity  of  pressure  implies  the  continuity  of  speed  across  it. 

Consider  a right-handed  system  of  rectangular  Cartesian  axes  with  Ox 
inclined  at  * small  angle  X to  the  undisturbed  stream  of  speed  U . Denote  the 
components  of  the  velocity  parallel  to  the  axes  Ox,  Oy  and  Oz  by 
U cos  X + u,  v and  w , so  that  u,  v and  w are  small  compared  with  U by 
the  assumptions  of  small  disturbances  and  small  inclination.  Suppose  the  vortex 
sheet  or  body  surface,  £ , is  defined  by  F(x,y,z)  - 0 and  also,  in  parametric 
form,  by  y - f (x,rO , z - g(x,ri)  , where  n is  a parameter  which  varies  along 
the  curve  W in  which  the  cross-flow  plane  x - constant  meets  £ . Then 

f(x,  f(x,n),  g(x,n))  & 0 , 


where  the  identity  sign  means  that  equality  holds  for  all  values  of  x and  n 
in  the  ranges  which  define  £ . Hence  the  identity  can  be  differentiated  with 
respect  to  x and  n to  give: 


F+Ff+Fg  -0 
x y x z°x 


F f + F g - 0 . 

y n z*n 


Since  £ is  a stream  surface  of  the  three-dimensional  flow: 


(A-1 ) 
(A-2) 


(U  cos  X + u)F  + vF  + wF  ■ 0 , 

x y z * 

or,  since  X is  small  and  u is  small  compared  with  U 
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UF„  + vF  + wF 
x y z 


0 . 


(A-3) 


The  three  equations  (A-l)  to  (A-3)  are  linear  and  homogeneous  in  F , F and 

w 

F , and  they  have  a non-trivial  solution,  so  the  determinant  of  the  coefficients  ; 

s 5 

vanishes.  Expanding  the  determinant  gives:  j 

3 


*£n  ' V8n  - °<Vx  * 8nfx> 


M 

i 


This  provides  an  equation  for  the  component  vq  of  the  velocity  in  the  cross* 

flew  plane,  normal  to  « . If  the  sense  of  the  normal  is  given  by  a positive  j 

(anticlockwise)  rotation  from  the  direction  of  increasing  n » j 

| 

5 

2 2 4 ‘i 

vn  - w cos  * - v sin  ♦ - (wf  - vg  )/(f  + g^)*  , ] 

| 

where  $ is  the  inclination  of  the  tangent,  in  the  sense  of  r\  increasing,  to  j 

\ 

the  axis  Oy  , and  where  the  positive  value  of  the  radical  is  understood.  Hence 


f g_  - g f 

nx  n x 

.2  . XT 


(C 


+ g ) 


(A-4) 


The  parameter  n does  not  appear  explicitly  in  (A-4) , but  the  same  parameter 
must  be  held  constant  in  the  two  partial  derivatives.  Particular  cases  of  (A-4) 
are 


3z 


COS  ijl 


iz 

3x 


sin  \|<  » 


3r  I . . 

' 3lle  *”  * 


(A-5) 


where  (r,6)  are  polar  coordinates  and  6 is  the  inclination  of  the  tangent  to 
the  radius  vector,  as  in  Fig  3. 

It  is  sometimes  convenient  to  work  with  the  stream  function,  ¥ , in  the 
cross-flow  plane.  If  the  velocity  potential  * is  the  real  part  of  the  complex 
potential  W(Z)  , where  Z ■ y + iz  , then  ¥ is  the  imaginary  part  of  W . 
Hence,  on  W 


3* 

3n 


4 


i 


3F 

Do 


■ - v 
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Hence,  by  (A-4) : 


Y - U J sin  Y - !“  cos  i}i)da  • U J dz  - |j  dy^  . (A-6) 


Equations  (A-4)  to  (A-6)  are  equivalent  expressions  of  the  condition  that  Z is 
a stream  surface. 

The  second  boundary  condition  applying  when  Z is  a vortex  sheet  is  that 
the  fluid  speed  is  continuous  across  Z , ie 


A cos  X + u)2  + v2  + w2^  - 0 , 


where  A is  a difference  operator  across  Z . To  be  definite,  A is  chosen  so 
that  AA  » Aj  - A^  , where  the  suffix  1 denotes  the  side  of  Z towards  which 
the  normal  points  and  2 the  opposite  side.  Because  X and  the  disturbances  are 
small,  this  condition  reduces  to 


2UAu  + A(v2  ♦ w2)  - 0 . 


(A-7) 


Now,  if  vfc  is  the  component  of  velocity  in  the  cross-flow  plane,  tangential  to 
, in  the  sense  of  n increasing. 


2.2 
v + w 


2 2 
Vt  + vn  5 


while,  to  the  order  of  accuracy  of  this  analysis,  (A-4)  shows  that  vq  is  con- 
tinuous across  Z . Hence  (A-7)  reduces  to 


UAu  + v Avfc  * 0 , 


(A-8) 


where  v is  the  mean  value  of  across  Z . Now,  if  the  velocity  potential 

tn 

in  the  fluid  on  the  two  sides  of  Z is  4>j(x,y,z)  and  » 


Au  * ^ 

x x 


*'y»* 


On  the  other  hand,  the  jump  in  $ across  Z is 


A*  - *i(x»  f(x»h)»  g(x,n))  “ *2 (x*  g(*»h)) 
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therefore 


A*  + f it  + g it 
x a y x s 

iu  + f iv  + g iw  . 
X ®x 


(A-9) 


Since  v is  continuous  across  Z , 
n 


Av  ■ Avfc  cos  ^ 


and 


Aw  * Avt  sin  ^ , 


so  that,  by  (A-9), 


Au 


3A« 


■gj-  - Avt(fx  cos  <|*  + gx  sin  *) 


Hence,  the  condition  of  continuity  of  pressure,  (A-8),  becomes 


9At  . I 3y  . ^ . . ml 

a r ■ 4vthi  <=o.  — 


(*-10) 


As  in  (A-5),  the  parameter  n does  not  appear  explicitly  in  the  boundary 
condition,  but  the  partial  derivatives  on  both  sides  of  (A-10)  must  be  calculated 
for  fixed  values  of  the  same  parameter.  If,  for  example,  the  polar  angle  0 is 
fixed,  (A-10)  becomes 


3A8 | . 3r I , Sal 

ir|9  ■ • 


(A-ll) 


For  conical  flow,  with  s • Kx  , this  reduces  to  the  familiar  form 

A*  ■ Av.  (x  cos  ♦ 


Similarly,  Clark's  equations,  (A-18)  and  (A-19)  of  Ref  5,  are  obtained  from 
(A-5)  and  (A-10)  when  the  sheet  is  defined  in  terms  of  his  parameters. 
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c,  d 
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u,  v,  w 


U 


*t  y»  * 

z 

z* 


LIST  OF  SYMBOLS 

diameter  of  circular  arc  in  Z^-plane 

undetermined  real  constants 

coordinates  of  points  B,  C,  D in  Z^-plane 

arbitary  positive  constant 

pressure  coefficient,  (p  - p„)/ipU^  • 

curve  in  which  £ meets  the  cross-flow  plane 

mapping  function 

parameter  describing  local  shape  of  wall,  see  (94) 
constants  defined  by  (14a) 
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pressure 
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time 

speed  of  undisturbed  stream 
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complex  potential 
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y + iz,  cross-flow  plane 
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6 

4 

c 

n 

e 

x 

e 

p 

o 

Z 

♦ 
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total  circulation  of  sheet 
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difference  operator  across  vortex  sheet 
small  complex  quantity,  see  (13) 

parameter  on  # in  Appendix;  dummy  variable  on  real  axis 
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inclination  of  x-axis  to  undisturbed  stream 
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velocity  potential 
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suffixes  S*  R denote  singular  and  regular  contributions 
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Fig  1 Section  of  wall,  separation  line  and  vortex  sheet,  showing  surface  streamlines 
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August  1967 

AD-907318.  Unclassified,  Dist.  E:  DoD  Only. 

Marking  pasted  on  document:  "Uncontrolled  copy.  Any  future  amendment 
will  be  issued  only  on  request  (DERA-DRIC  Kentigern  House  Glasgow). 

Again  DTIC  has  this  document.  What  does  "Uncontrolled  copy"  mean;  does 
it  change  the  distribution  audience? 

3)  This  may  be  for  Chrissie  Stewart's  attention.  Let  me  know  if  you  want 
me  to  fax  copies  of  the  dueling  RDPs 

DERA/CIS/CIS1/CR990854 

Evaluation  of  Modern  HF  Ray  Tracing  (Final  Report) 

Neil  Rogers 
November  1999 

The  U.S.  European  Office  of  Aerospace  Research  & Development(EOARD) 
http://www.ehis.navy.mil/eoardwebsite2/  sent  a copy  of  the  above  contracted 
report.  EOARD  completed  a RDP  which  states  that  the  report  is  "Public 
Release,  Distribution  Unlimited."  However,  DERA  marked  every  page  "UK 
Unclassified  - US  for  Official  Use  Only."  EOARD  is  adamant  that  it  and  not 
DERA  is  the  controlling  office  and  that  DERA  does  not  have  the  right  or 
authority  to  mark  the  document  on  behalf  of  the  U.S.  Air  Force. 

EOARD  should  probably  have  handled  this  by  asking  DERA  to  produce  a copy 
of  the  report  without  the  markings.  However,  they  didn't,  and  I have  the 
document  here  waiting  to  be  processed. 

Can  I get  an  email  or  statement  from  DRIC  or  DERA  stating  that  this 
document  done  under  U.S.  Air  Force  Contract  F61775-98-W-WE133  is  Public 
Release  as  determined  by  EOARD? 

Bonnie  Klein 

Program  Manager  for  Foreign  Acquistions 
Defense  Technical  Information  Center 
Ft.  Belvoir,  VA 
(703)  767-8037 


The  Information  contained  in  this  E-Mail  and  any  subsequent  correspondence 
is  private  and  is  intended  solely  for  the  intended  recipient(s). 

For  those  other  than  the  recipient  any  disclosure,  copying,  distribution, 
or  any  action  taken  or  omitted  to  be  taken  in  reliance  on  such  information  is 
prohibited  and  may  be  unlawful. 
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